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Abstract. In this paper we study, from a variational and geometrical point 
of view, second-order variational problems on Lie groupoids and the construc¬ 
tion of variational integrators for optimal control problems. First, we develop 
variational techniques for second-order variational problems on Lie groupoids 
and their applications to the construction of variational integrators for opti¬ 
mal control problems of mechanical systems. Next, we show how Lagrangian 
submanifolds of a symplectic groupoid gives intrinsically the discrete dynamics 
for second-order systems, both unconstrained and constrained, and we study 
the geometric properties of the implicit flow which defines the dynamics in the 
Lagrangian submanifold. We also study the theory of reduction by symmetries 
and the corresponding Noether theorem. 


1. Introduction. The topic of discrete Lagrangian mechanics concerns the study 
of certain discrete dynamical systems on manifolds. As the name suggests, these 
discrete systems exhibit many geometric features which are analogous to those in 
continuous Lagrangian mechanics. In particular, the discrete dynamics are de¬ 
rived from variational principles, have symplectic or Poisson flow maps, conserve 
momentum maps associated to Noether-type symmetries, and admit a theory of 
reduction. While discrete Lagrangian systems are quite mathematically interest¬ 
ing, in their own right, they also have important applications in the design of 
structure-preserving numerical methods for many dynamical systems in mechanics 
and optimal control theory. 

Numerical methods which are constructed in this way are called variational in¬ 
tegrators. This approach to discretizing Lagrangian systems was put forward in 
papers by Bobenko and Suris [6], Moser and Veselov [60], and others in the early 
1990s, and the general theory was developed over the subsequent decade (see Mars- 
den and West [51] for a comprehensive overview). 

2010 Mathematics Subject Classification. Primary: 70G45 ; Secondary: 70Hxx, 49J15, 53D17, 
37M15. 

Key words and phrases. Discrete Lagrangian mechanics, higher-order variational problems, 
optimal control, Lagrangian submanifolds. Lie groupoids, Lie algebroids, geometric integration. 


1 



2 


LEONARDO COLOMBO AND DAVID MARTIN DE DIEGO 


A. Weinstein [62] observed that these systems could be understood as a special 
case of a more general theory, describing discrete Lagrangian mechanics on arbitrary 
Lie groupoids where the Lagrangian function is defined on a Lie groupoid. A Lie 
groupoid G is a natural generalization of the concept of a Lie group, where now 
not all elements are composable. The product gh of two elements is only defined 
on the set of composable pairs G 2 = {{g^h) S G x G\[i[g) = a{h)}, where a : 
G ^ Q and (5 : G ^ Q are the source and target maps over a base manifold 
Q. This concept was introduced in differential geometry by Ereshmann in the 
1950’s. The infinitesimal version of a Lie groupoid G is the Lie algebroid tag ■ 
AG Q, which is the restriction of the vertical bundle of a to the submanifold of 
the identities. This setting is general enough to include the discrete counterparts of 
several types of fundamental equations in Mechanics as for instance, Euler-Lagrange 
equations for Lagrangians defined on tangent bundles [51], Euler-Poincare equations 
for Lagrangians defined on Lie algebras [49], [50], Lagrange-Poincare equations for 
Lagrangians defined on Atiyah bundles, etc. Such discrete counterparts are obtained 
discretizing the continuous Lagrangian to the corresponding Lie groupoid and then 
applying a suitable discrete variational principle. 

A complete description of the discrete Lagrangian and Hamiltonian mechanics 
on Lie groupoids was given in the work by Marrero, Martin de Diego and Martinez 
[41] (see also [42]). Following the program proposed by A. Weinstein [62], in this 
work, we generalize the theory of discrete second-order Lagrangian mechanics and 
variational integrators for a second-order discrete Lagrangian : G 2 —> K in two 
main directions. First, we develop variational principles for second-order variational 
problems on Lie groupoids and we show how to apply this theory to the construction 
of variational integrators for optimal control problems of mechanical systems. Sec¬ 
ondly, we show that Lagrangian submanifolds of a symplectic groupoid (cotangent 
groupoid) give rise to discrete dynamical second-order systems. We also develop 
a reduction by symmetries, and study the relationship between the dynamics and 
variational principles for second-order variational problems. Finally we study dis¬ 
crete second-order constrained Lagrangian mechanics on Lie groupoids. The main 
application of this theory will be systems subjected to constraints and underactu¬ 
ated control systems. Our results in the second part of the paper are based on the 
papers [41] and [43] but in second-order theories. 

There are variational principles which involves higher-order derivatives [4], [9], 
[22] , [23] , [24] , [25] , [39] , [47] , [52] , [54] , [55] since from it one can obtain the equations 
of motion for Lagrangians where the configuration space is a higher-order tangent 
bundle. The study of higher-order variational systems has regularly attracted a lot 
of attention from the applied and theoretical points of view (see [37]), but recently, 
higher-order variational problems have been studied for their important applications 
in optimal control, aeronautics, robotics, computer-aided design, air traffic control, 
trajectory planning and computational anatomy. 

Since the main applications of higher-order variational problems are second-order 
problems we will focus our attention in the second-order case along the work, leaving 
the extension to higher-order case as a straightforward development. 

The organization of the paper is as follows. In Section 2 we recall some construc¬ 
tions and results on discrete Mechanics on Lie groupoids which will be used in the 
next sections. Section 3 is devoted to study variational principles for second-order 
discrete mechanical systems on Lie groupoids, their extension to the constrained 
case and the application to the theory of optimal control of mechanical systems 
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and construction of variational integrators. In Section 4 we show how Lagrangian 
submanifolds of an appropriate symplectic groupoid (cotangent groupoid) give rise 
to discrete dynamical second-order systems. From such Lagrangian submanifold 
we obtain the discrete second-order Euler-Lagrange equations on Lie groupoids and 
such equations correspond with the ones obtained from the variational point of view. 
We also develop a reduction by symmetries, and we study the relationship between 
the different dynamics and variational principles for these second-order variational 
problems. Finally we study discrete constrained second-order Lagrangian mechan¬ 
ics. This allows for systems with arbitrary constraints. 

Throughout the paper, we have occasion to draw on certain technical construc¬ 
tions using the theory of Lie algebroids and retraction maps. We provide some 
supplementary details and discussion of these and an overview on discrete mechan¬ 
ics and higher-order tangent bundles in some Appendices at the end of this paper. 

2. Groupoids and discrete mechanics. This section review some results about 
Lie groupoids and discrete mechanics on Lie groupoids based on [41] and [42]. 

2.1. Generalities about Lie groupoids. A groupoid is a small category in which 
every morphism is an isomorphism (i.e. all morphism is invertible). That is. 

Definition 2.1. A groupoid G over a set Q, denoted G ^ Q, consists of a set of 
objects Q, a set of morphisms G, and the following structural maps: 

• a source map a: G ^ Q and a target map /3: G ^ Q. Thus an element g € G 
is thought as an arrow from x = a{g) to y = l3{g) in Q. 



X = a(g) y = /3{g) 

• an associative multiplication map m: G2 —>■ G, m{g,h) = gh, with a{g) = 
a{gh) and /3(/i) = (i{gh) where 

G2-.= Gpx^G = {{g, h)GGxG\ p{g) = a{h)} 

is called set of composable pairs defined by a and /3. gh is thought as the 
composite arrow from a; to z if 5 is an arrow from x = ot{g) to y = Pig) = a(/i) 
and h is an arrow from y to /3(h) = 0 . 



X = a(g) = a(gh) y = p{g) = a(h) z = g(h) - P{gh) 

• an identity map, e: Q —)■ G, a section of a and /3, such that for all g G G, 

e(a(5))5 = 5 = 5^ iPig)), 

• an inversion map i: G ^ G, mapping g into g~^, such that for all g G G, 

gg-^ = e (a(g)) and g-^g = e (P(g)). 



X = a(g) = p(g L y = P{g) = a(g L 



Remark 1. Alternatively, a groupoid can be seen as a weak version of a group, 
where the multiplication will be defined only for elements in G 2 C G x G. o 
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We will focus on a particular class of groupoids, the Lie groupoids which have a 
differential structure in addition to their algebraic structure. 

Definition 2 .2. A Lie groupoid is a groupoid G ^ Q where 

1 . G and Q are differentiable manifolds, 

2 . a,/3 are submersions, 

3. the multiplication map m, inversion i, and identity e, are differentiable. 

Remark 2. In Definition 2.2, a and (3 must be submersions so that G 2 is a dif¬ 
ferentiable manifold. From the definition it follows that m is a submersion, e is an 
immersion, and I is a diffeomorphism. o 

If g e Q, oi~^{q) and will be said the a-fiber and ,5-fiber of q. 

Definition 2.3. Given a groupoid G ^ Q and g G G, define the left translation 
ig-. a~^ (/3(g)) —t a.~^ (o-(g)) and right translation rg-. j3~^ ( 0 ( 5 )) —>■ f3~^ (/^(s)) by 
g to be 

ig(h) = gh and rg(h) = hg. 

Note that, = ig-i and r~^ = r^-i. 

Denoting by X(G) the set of vector fields on G one may introduce the notion of 
a left (right)-invariant vector field in a Lie groupoid, as in the case of Lie groups. 

Definition 2.4. Given a Lie groupoid G Q, a vector field X G X(G) is left- 
invariant (resp., right-invariant) if X is a-vertical (resp., /3-vertical), that is, it 
is tangent to the fibers of a (resp., /3), Ta(X) = 0 (resp., Tf3(X) = 0) and 
(Th(.g) (X(h)) = X (gh) for all (g,h) G G 2 (resp., (Tnrg) (X(h)) = X (hg)). 

Definition 2.5. A Lie algehroid A over a manifold Q is a real vector bundle ta '■ 
A —>• Q together with a Lie bracket I-, •] on r(r^), the set of sections of : A —>• Q, 
and a bundle map p : A —>• TQ such that [AT,/F] = /|Ar, F] -|- p(X)(f)Y for all 
X,Y G r(r^) and / S G°°(Q). 

Remark 3. If (A, l-,-],^) is a Lie algebroid over Q then p is a homomorphism 
between the Lie algebras (r(ryi), I-, •]) and (X(Q), [•, •]). o 

In Lie groups, the infinitesimal version of a Lie group is a Lie algebra, therefore 
we will see that the corresponding infinitesimal version of a Lie groupoid is a Lie 
algebroid. Next, we define the Lie algebroid associated with a Lie groupoid G ^ Q. 
Given a Lie groupoid G ^ Q, consider the vector bundle 

tag ■ AG — Q 

whose fiber at a point q G Q is AgG = kerTj(q)a = Fa, i.e., the tangent space to 
the a-fiber at the identity section, for q G Q. 

It is easy to prove that there exists a bijection between the space of sections 
r (tag) and the set of left-invariant vector fields on G. If A" is a section of tag, the 
corresponding left-invariant vector field on G will be denoted by ^ G X(G) where 

^(g)=(T,(P(g))£g)(X(/3(g))). (I) 

The Lie algebroid structure on AG is given by the bracket I-, •] and the anchor 
map p defined as follows: 

= [^, , p(X)(q) = (r,(,)/3) (X(q)) , 

for all X,Y gT (tag) and q G Q. 


( 2 ) 
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Definition 2.6. Given a Lie groupoid G ^ Q, the triple {AG, I-, -IjP) defined in 
(2) is called Lie algebroid associated to G ^ Q. 

Remark 4. Alternatively one can also establish a bijection between the space of 
sections L {tag) and the set of right-invariant vector fields on Q, by 

^( 5 ) = - {Te{aig)) {rg ° i)) (a(5))), (3) 

which yields the Lie bracket relation 

Thus the mapping X into is a Lie algebra isomorphism, and the mapping X into 
is a Lie algebra anti-isomorphism (see [20] and [40] fore more details). o 

2.1.1. Examples of Lie groupoids: We introduce some examples of Lie groupoids. 

• The pair or banal groupoid: Let Q be a differentiable manifold, and consider 
the product manifold G = Q x Q. Gisa Lie groupoid over Q where the source 
and target maps a and /3 are the projections onto the first and second factors 
respectively. The identity is defined as e{q) = {q, q) for all q G Q, the multiplication 
m{{q, s), {s, r)) = {q, r) for {q, s),{s,r) G Q x Q and the inverse map i{q, s) = {s, q). 

Note that, if g is a point of Q, then ~ TqQ. Hence the Lie algebroid of 

G is isomorphic to the standard Lie algebroid ttq : TQ Q. In this sense, the 
Banal groupoid is considered as the discrete space for discretizations of Lagrangian 
functions L : TQ ^ M.. 

• Lie groups: Let G be a Lie group. G is a Lie groupoid over one point Q = {e}, 
the identity element of G. The structural maps of the Lie groupoid G are 

“( 5 ) = e. I3{g) = e, e(e) = e, i{g) = g~^ and m{g, h) = gh, for g,hGG. 

The Lie algebroid associated with G is the Lie algebra g of G. 

• Transformation or action Lie groupoid. Let H he & Lie group with identity e 
and tf ■. Q X H ^ Q he a. right action of on Q- The product manifold G = Q x H 
is a Lie groupoid over Q, with structural maps given by 


a{q,h)=q, f3{q,h) = ip{q,h), e{q) = {q,e), m{{q,h), {(p{q,h),h')) = {q,hh'), 

i{q, h) = {(p{q, h),h~^) ioi q G Q and h, h' G H. 

The Lie groupoid G is called action or transformation Lie groupoid and its associated 
Lie algebroid is the action algebroid pr^^ : M x t) —>■ M where f) is the Lie algebra of 
the Lie group H (for more details, see [40] and [41]). 

• The cotangent groupoid: Let G Q be a Lie groupoid. If A*G is the dual 
bundle to AG then the cotangent bundle T*G is a Lie groupoid over A*G. The 
projections j3 and a, the partial multiplication ©t*Gj the identity section e and the 
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inversion i are defined by the structural maps of G as follows, 

= Ms(m(/ 3 ( 3 )) 4 )(^)). for e T;G and X G 
dK)(y) = (r,(„(;,))(eo/3))(y))), 

for Vh^T^G and Y G Ao,(h)G, 

{Hg ©t*g ’^h){T(g^h)'m{Xg, Yh)) = tJ.g{Xg) + Vh{Yh), 

for {Xg,Yh)GT(^g,h)G 2 . 

l^x{,Xf^(^x) (^e(a:) (c ® Qi)) ))), 

for € AxG and X^f^x) ^ Y^(^x)G^ 

\{^Xg){Xg-^) = -^Jig{{Tg-^^){Xg-^)), for G T*gG and X^-i G T^-iG. 

Here a,/3,m,f,e are the structural maps of G =| Q (for more details, see [20] and 

[31]). 

• Symplectic groupoids: Finally, we introduce a subclass of Lie groupoids with 
an additional structure, symplectic groupoids. They are endowed with a symplec¬ 
tic manifold structure. A symplectic groupoid is a Lie groupoid G ^ Q with a 
symplectic form w on G such that the graph of the composition law m given by 

graph(m) := {(g, /i, r) G G x G x G | {g, h) G G 2 and r = m{g, h)} 

is a Lagrangian submanifold of G x G x G~ with the product symplectic form, where 
the first two factors G are endowed with the symplectic form ui and the third factor 
G~ with the symplectic form —oj. 

Observe that if G =t Q is a Lie groupoid then the cotangent groupoid T*G ^ 
A*G is a symplectic groupoid with the canonical symplectic 2-form on T*G, denoted 
log- 

2.2. Lie Groupoids and Discrete Mechanics. Next, we give a review of some 
generalities on discrete mechanics on Lie groupoids based in [41] and [42] . 


2.2.1. Discrete Euler-Lagrange equations. Let G be a Lie groupoid over Q with 
structural maps 


a, (3 : G —>■ Q, e : Q —G, i : G —G, m : G 2 


G. 


Denote by tag ■ AG —>■ Q the Lie algebroid of G. 

A discrete Lagrangian is a function : G —>■ M. Fixed g G G, we define the set 

of admissible sequences with values in G: 

Cf = G G^ I ( 5 fc, 5 ffc+i) G G 2 for fc = 1,... ,X - 1 and 51 .. . 57 V = g}. 

(5) 

An admissible sequence {gi,..., gw) G is a solution of the discrete Euler- 
Lagrange equations if 


0 = 


IV-l 

k{gk){^d) — ^k{gk-ei)(^d) , 


k=l 


for Xi, 


• jXjv-i G r(r_AG). 


For N = 2 we obtain that ( 5 , h) G G 2 is a solution of the discrete Euler-Lagrange 
equations if 


y(5)(Ld)-y(/i)(Ld) = 0 


for every section X of AG. 
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Remark 5. Marrero et al. [41] showed that these discrete Euler-Lagrange equations 
are also equivalent to the sequence gi,... G G corresponding to a critical point 
of the action sum 

N 

{gi,...,gN) I—^ ^Ld(gfc), 

fc=i 

over the space of admissible sequences. 

In the case when G is the banal groupoid Q x Q ^ Q, this recovers the discrete 
Euler-Lagrange equations, 

Dlhdi<lk, Qk+i) + D2Ld{qk-i,qk) = 0 

for fc = 1,..., TV — 1, as in Marsden and West [51]. o 

2.2.2. Discrete Lagrangian evolution operator. We say that a differentiable mapping 
Ip : C? —>■ G is a discrete flow or a discrete Lagrangian evolution operator for if 
it verifies the following properties: 

- graph(^') C G 2 , that is, (g, ^'( 5 )) e G 2 , V 5 S G. 

- ( 5 , 41 ( 5 )) is a solution of the discrete Euler-Lagrange equations, for all g G G, 
that is, 

^(g)(Ld)-^(vI/( 5 ))(Ld) = 0 ( 6 ) 

for every section X of AG and every g G G. 

2.2.3. Discrete Legendre transformations. Given a discrete Lagrangian : G —>■ K 
we define two discrete Legendre transformations F“Ld : G -G A*G and F+L^ : G -G 
A*G as follows (see [41]) 

(F"Ld)(/i)(ne(a(/i))) = - v^(^a(h))(.^d orhO i), for Pe(a(/i)) e Aa(h)G, (7) 
(F l^d){9){'^e{/3{g))) O £g), for G Ap(^g'^G. ( 8 ) 

Note that (F+Ld)(g) G ^p(g)G and (F“Ld)(Ii) G A*^f^f^^G. 

2.2.4. Regular discrete Lagrangians and Hamiltonian evolution operator. A discrete 
Lagrangian : G —>■ K is said to be regular if and only if the Legendre transfor¬ 
mation F“'"Ld is a local diffeomorphism (equivalently, if and only if the Legendre 
transformation F“Ld is a local diffeomorphism). In this case, if {go, Hq) G G x G is 
a solution of the discrete Euler-Lagrange equations for then, one may prove (see 
[41]) that there exist two open subsets Uq and Vq of G, with go G Uo and ho G Vo, 
and there exists a (local) discrete Lagrangian evolution operator 41 : Go —>■ Vq such 
that: 

1 - ^( 50 ) = ho, 

2. 41 is a diffeomorphism and 

3. 4> is unique, that is, if Uq is an open subset of G, with go G Uq, and 41 ' : Gq —>■ G 
is a (local) discrete Lagrangian evolution operator then 

^ lc/onc/'= • 

Moreover, if F+L^ and F“Ld are global diffeomorphisms (that is, is hyperregular) 
then 41 = (F“Ld)“^ oF+L^. 

If : G —> K is a hyperregular Lagrangian function, then pushing forward to 
A*G with the discrete Legendre transformations, we obtain the discrete Hamiltonian 
evolution operator, 41 ■. A* G ^ A* G given by 

$ = F±Ld o 4- o (F±Ld)-^ = F+Ld o (F-L^)"!. 


(9) 
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2.2.5. Example: Discrete Euler-Poincare equations. Let G be a Lie groupoid over 
Q = {e}, the identity of G. Given ^ G g we have the left and right invariant vector 
fields 




e ( 5 ) = (re4)(C) and ^ (g) = (T,rg)(C), for g G G. 

Given a discrete Lagrangian : G —)■ M its discrete Euler-Poincare equations 


are 


(T’e4J(e(Ld) - (T’erg^+J(0(Ld) = 0, for all ^ e 0 and gk,gk+i £ G, 
that is, 

(^;,dLd)(e) = (rgV^dL,)(e) 

(see [6, 49, 50]). 

A way to discretize a continuous problem is by using a retraction map r : g —> G, 
which is an analytic local diffeomorphism and maps a neighborhood 14 C g of 0 £ g 
to a neighborhood G C G of the identity e £ G. We have that = e for all 

C £ g (see [7]). 

The retraction map provides a local chart on the Lie group and it is used to 
express a small change in the group configuration through a unique Lie algebra 
element, namely where > 0 is a small enough time step, 

Cfc £ 0j £ r C g and gk,gk+i £ G C G, i.e., if ^k were regarded as an average 
velocity between gk and gk-^-l^ then r is an approximation of the corresponding 
vector field on G (see [33]). 

To derive the discrete Euler-Poincare equations, one uses the left-trivialized tan¬ 
gent retraction map dr^ : 0 —t g and its inverse drf^ : 0 —t g defined by 

T^T{r]) =T,irioi^^dv)) = dT5(7?)T(C), 

for 77 G g (see [32], [7]) . 

The Lie algebra g is on itself a vector space, then it is natural to consider local 
coordinates on g. We will write T{hf) = g for a enough small time step h > 0 such 
that C V where E is a local neighborhood of 0 G g. Fixing a basis {e-y} of g 
we induce coordinates (j/'*') on g. In these coordinates, a basis of left-invariant and 
right-invariant vector fields is 

= Tr(hn)T~^{T{hvi)e-^) = dTff.^{Adr(hri)e-^), 

= Tr(hr,)r~He'fT{hr])) = dr-^^ej), 

where ?? G g. 

Given a Lagrangian G g —> K, the discrete Euler-Poincare equations are: 

KiVk)il)-'^iVk+i){l) = 0, ( 10 ) 

that is, 

(see [42]). 
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3. Second-order variational problems on Lie groupoids and optimal con¬ 
trol applications. In this section, we discuss discrete second-order Lagrangian 
mechanics using techniques of variational calculus on Lie groupoids (see [30] and 
[41] for first order variational calculus on Lie groupoids) and we illustrate our results 
with some examples and applications in the theory of optimal control of mechanical 
systems. 

3.1. Second-order variational problems on Lie gronpoids. Let G be a Lie 

groupoid with structural maps a, /S : G ^ Q, e : Q ^ G, i : G ^ G and m : G 2 —>■ 
G. Denote by tag ■ 4lG —>■ Q the Lie algebroid associated with the Lie groupoid G. 

Definition 3.1. A discrete seeond-order Lagrangian Ld '■ G 2 —>■ K is a differentiable 
function defined on the set of composable elements describing the dynamics of the 
mechanical system. 

We denote by G^ the product G x G x G x G. As in the first order case, fixed 
g G G, we define the set of admissible sequences in G^ with values in G by consider 
A^ = 4 in (5), that is, 

Cg = {( 51 . 52 , 93 , 94 ) e G^ I {gk,gk+i) G G 2 for A: = 1,2,3 
with gi and 54 fixed and 51329354 = 51- 

Given a tangent vector at the point 5 = ( 51 , 92 , 93 , 54 ) to the manifold Gg, we 
may write it as the tangent vector at t = 0 of a curve c{t) in G^, which passes 
through 5 at t = 0. This type of curves has the form 

c{t) = (5i,92/i2(A),^^^(A)93,54), (11) 

where /i 2 (t) G a~^{P{g 2 )), for all t, and /i2(0) = e{(3{g2))- The curve c is called a 
variation of 5 . Therefore we may identify the tangent space to Gg at 5 with 

TgGg = {v 2 I V 2 G A^^G where X 2 = I3{g2)}- 

The curve V 2 is called infinitesimal variation of 5 and is the tangent vector to the 
a-vertical curve /12 at t = 0. 

We define the discrete action sum associated with the discrete second-order La¬ 
grangian : G 2 —>■ K by 

K 
3 

^ , ^djgkj gk+i)- (12) 

fc=i 

To derive the discrete equations of motion we apply Hamilton’s principle of crit¬ 
ical action. In order to do that, we need to consider the variations of the discrete 
action sum. 

Definition 3.2 (Discrete Hamilton’s principle on Lie groupoids). Given 5 S G, an 
admissible sequence 5 G Cg is a solution of the Lagrangian system determined by 
Ld : G 2 —>■ K if and only if 5 is a critical point of Sl^ ■ 

Proposition 1. Given g G G, the admissible sequence 5 = ( 51 , 52 , 53 , 54 ) G Gg is 
a solution of the Lagrangian system determined 65 Ld : G 2 —>■ M if and only if 5 
satisfies the discrete second-order Euler-Lagrange equations for Ld : G 2 —>■ K given 
by 

tg^ (L>iLd( 52 , 53 ) + L) 2 Ld( 5 i, 52 )) + {rg^ oi)* (DiLd( 53 , 54 ) + L) 2 Ld( 52 ,53)) = 0 . 

(13) 


Sl, : G" 
9 
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Proof. By definition ( 3 . 2 ), 5 is a solution of the Lagrangian system determined 
by Ld : G2 —> K if it is a critical point of Sl^- In order to characterize the critical 
points, we calculate. 


d 

dt 


SrAcit)) = -r, 

t=o dt 


t^o 


{Ld{gi,g2h2it))+Ld{g2h2{t),h2 {t)g3)+Ldih.^ (^ 53 , 54 )} 


where c(<) is the variation of g defined in (11). 

Then, the condition — 

’ dt 


t=o 


SL^(c(t)) = 0 is equivalent to 


0 =d(Ld o 4J(e(/3(5f2)))(i'2) + d(Ld o rg^ o i)(e(/3(g2)))(f2) 

+ d{Ld o Tg^ o z)(e(/3(53)))(u2) + d{Ld o eg2){e{/3{gi))){v2), 


( 14 ) 


where V2 G is the infinitesimal variation of g, e{j 3 {g 2 )) = li2(0) and £g and 

Tg were defined on Definition 2 . 3 . 

Therefore, 5 is a solution of the Lagrangian system determined by the discrete 
second-order Lagrangian : G2 —>■ K if and only if it satisfies equations ( 14 ), that 
is, 5 is a solution of the Lagrangian system determined by : G2 —>■ K if and only 
if g satisfies 


^92 iDiLd{g 2 ,g 3 ) + D 2 Ld{gi,g 2 )) + {rg^ o i)* {DiLd{g 3 ,g 4 ) + D 2 Ld{g 2 ,g 3 )) = 0. 

( 15 ) 

□ 

The equations given above are called discrete second-order Euler-Lagrange equa¬ 
tions. 


3 . 1 . 1 . Example: Discrete second-order Euler-Lagrange equations on the pair groupoid. 
Let Q X Q =$ Q be the Banal groupoid. An admissible path is the 4 -tuple 
((91,92), (92,93), (93,94), (94,®)) e (Q ^ <3)2 is isomorphic to 3 Q where 

inclusion of 3 Q into (Q x (5)2 is given by the map (q, q, q) ^ ((q, q), (q, q)). Apply¬ 
ing Hamilton’s principle for the discrete second-order Lagrangian Ld : (Q x (5)2 — 
3 Q -)-R given by Ld((qfe, qfc+i), (qfc+i, 9fc-i-2)) = Ld(qfe, qfc+i, 9fe-i-2), one gets 


d 

dt 


t=o 


Sl,= 


3 

V 

‘=°^1 


.bd(qfe, qfc+1, 9 ^+ 2 ). 


Hence, the path ((91,92), (92,93), (93,94), (94,95)) G is a critical point of Sl^ 

if and only if it satisfies 


- 03 ^( 1 ( 91 , 92 , 93) + A' 2 Td( 92 ,q 3 , 94 ) + T>iTd( 93 , 94 , 95) = 0 . ( 16 ) 

for (91,92) and (94,95) fixed points in the Banal groupoid. 

These equations are the discrete second-order Euler-Lagrange equations for Ld : 
Q X Q X Q ^ R (see for example [ 1 ]). 


3 . 1 . 2 . Example: Discrete second-order Euler-Poincare equations. Let G be a Lie 
group, that is G is a Lie groupoid over the identity element {c} of G. Given 
9 G G and a discrete second-order Lagrangian Ld : G2 ^ K, the solution for the 
Lagrangian system determined by the discrete Lagrangian Ld are 

0 =^*gk^i^d{gk,gk+i) + (*g^D2Ld{gk-i,gk) - r*g^^^D2Ld{gk,gk+i) 

~ f*gk+i^^^d^9k+l, gk+ 2 ) 

for an admisible sequence {gk-i,gk,gk+i,gk+2) G with g = gk-igkgk+igk+2 and 
91,94 fixed points in G. 







SECOND-ORDER VARIATIONAL PROBLEMS ON LIE GROUPOIDS 


11 


The equations given above are the discrete second-order Euler-Poincare equations 
(see for example [17] and [9]). 

3.1.3. Example: Discrete second-order Euler-Lagrange equations on an action Lie 
groupoid. Let H he & Lie group, and Q a differentiable manifold. Let tp : QxH ^ Q 
be a right action, ip{q, h) = qh. We consider the action Lie groupoid, G = Q x H 
over Q. The set of composable elements is determined by 

G 2 = {{{qi,hi), (g 2 , / 12 )) I q 2 = qihi}. 

If {q,h) S G, the left-translation (-(q^h) ■ Oi~^{qh) —>■ a~^{q) and the right- 
translation /i) : j3~^{q) —)■ P~^{qh) (where a and jd are the source and target map 
of G) are given by £(,_,j)(g/i, h') = {q,hh') and r(q^h){q{h')-'^,h') = {q{h')-^,h'h) 
for q G Q,h,h' G H. 

Consider an admissible path 

[{qi,hi), (gi/ii, / 12 ), {qihih 2 ,h‘i), (gi/ii/ 12 / 13 , ^ 4 )) € C'^ 

where x = (qi, h) G Q x H with h = and {qi, hi), ((lihih 2 h^, ft. 4 ) are fixed 

in (Q X iJ). A discrete second-order Lagrangian is defined as L^ : {Q x H )2 — 
QxiJxiL—T^Kby 

Ldihi, hi, /i2) := Ldiiqi, hi), [qihi, /i2)). 

The discrete Euler-Lagrange equations for the system determined by the discrete 
second-order discrete Lagrangian Ld '■ Q x El x H ^ are determined by 


0 —G(^q-^hiM) {h)iLd{qihi, h2, /13) -b D2Ld{qihi, /12, ^3) + D^Ld{qi, hi, /12)) 

+ i'^(qihih2,h3) ° *)* {DiLd{qihih2, /13, ^4) + D2Ld{qihih2, /13, ^4) 
pD^Ldiqihi, h2, /13)) . 

92 =(p{qi,hi). 

3.2. Second-order constrained variational problems on Lie groupoids. 

Next, we extend the previous variational principle to second-order variational prob¬ 
lems for systems subject to second-order constraints. The constructions presented 
here are interesting for applications in optimal control problem of underactuated 
mechanical controlled systems. 

Let Ld : G 2 —>■ K be a discrete second-order Lagrangian describing the dynamics 
of a discrete mechanical system. Suppose that the dynamics is restricted. This 
restriction is given by the vanishing of s smooth constraint functions : G 2 —>■ K, 
A = 1,.. ., s determining a submanifold M of G 2 . 

The dynamics of the second-order constrained variational problem associated 
with Ld and is described by the discrete constrained second-order Euler-Lagrange 
equations determined by considering the augmented Lagrangian : G 2 x K® —>• K 
given by Ld = Ld -\- ^A^d where A^i = (Ai,..., As) G K® are Lagrange multipliers 
to be determined (see subsection 4.5 for an intrinsic approach). 

Given g G G, the set of admissible sequences is given by 

Gg,M = {(fli: 92 ,93j 94 , A^, A^, A^) S G"^ X K^® | {gk,gk+i) G G2,^d{9k,gk+i) = 0 

for k = 1,2,3 with gi and gi fixed and 91329394 = 9 }- 
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Consider the extended action sum associated with the extended Lagrangian : 
G 2 X K® -)■ K 

3 

( 5052 , 53 , 54 , I —!• ^ [Ld(5fc,5fe+i) + (A^O^'J’d (5fe,5fe+i)] ,(17) 

fc=i 

where A^ = (Af,...,Aj) e K®. An easy adaptation of the variational principle 
(3.2) for the discrete extended Lagrangian Ld can be done to obtain the discrete 
constrained second-order Euler-Lagrange equations by extremizing the extended 
action sum . The equations describing the dynamics of second-order constrained 
variational problems are 

( 51 , 52 ), 0 = $;f(g2,53), 0 = $;^( 53 ,g 4 ) with A = 1 ,... ,s, 

*1 =^92 (Tilled(52,53) + A^I?id);^(g2,53) + D2Ld{gi,g2) + A)iD2‘hd (51,52)) ( 18 ) 

+ (^93 ° *)* {DiLdigs,gi) + A^£>i$d (53,54) + D2Ld{g2, gs) + A^Ll2‘l’d (52,53)) ■ 

3.3. Application to optimal control of mechanical systems. In this section 
we study how to apply the second-order Euler-Lagrange equations on Lie groupoids 
to optimal control problems of mechanical systems defined on Lie algebroids. After 
introducing optimal control control problems, we study their discretization. 

3.3.1. Optimal control problems of total-actuated mechanical systems on Lie alge¬ 
broids. Let {A, I’, ’I, p) be a Lie algebroid over Q with bundle projection ta ■. A ^ Q. 
The dynamics is specified fixing a Lagrangian L : A —)■ M (see Appendix C) . Ex¬ 
ternal forces are modeled, in this case, by curves up : K —t A* where A* is the dual 
bundle : A* —>■ Q. 

Given local coordinates ( 9 ®) on Q, and fixing a basis of sections {ca} of : A —>• 
Q we can induce local coordinates ( 5 *, 5 “) on A; that is, every element b G Aq = 
(q) is expressed univocally as b = y°'ea{q). The notion of admissible curves 
replaces that of natural prolongation in the context of Lie algebroids. 

Definition 3.3. Let (A, I-, •], p) be a Lie algebroid over Q with projection ta : A —>■ 
Q. A curve ^ / C K —t A is an admissible curve on A if 

In a local description, a curve f on A given by fit) = {q"^ (t), y°‘(t)), is admissible 
if 

5* = 5L(5)y“ 

d 

where if & = y^Caiq) with q = tag(&) then p{b) = Paiq)y°‘^ ■ 

oq^ 9 

It is possible to adapt the derivation of the Lagrange-d’Alembert principle to 
study fully-actuated mechanical controlled systems on Lie algebroids (see [19] and 
[48]). Let qQ and qp fixed in Q, consider an admissible curve f : I cM. —>■ A which 
satisfies the principle 

0 = 6 f L{f{t))dt-\- ( {uFit),f]{t))dt, 

Jo Jo 

where p{t) G A.r^(j(t)) and upit) £ ^T^(^(i)) defines the control force (where we 
are assuming they are arbitrary). The infinitesimal variations in the variational 
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principle are given by = rf^ ^ for all time-dependent sections r] G r(ryi), with 
77 ( 0 ) = 0 and r)(T) = 0, where is a time-dependent vector field on A, the 
complete lift, locally defined by 

(see [19], [44], [45] and [46]). Here the structure functions are determined by 
|e/3) e-yj Cp^Ca. 

From the Lagrange-d’Alembert principle one easily derives the controlled Euler- 
Lagrange equations by using standard variational calculus 

d f dL\ ^ dL a dL 


The control force mf is chosen such that it minimizes the cost functional 


[ C{q\y°‘,{uF)a)dt, 

Jo 

where C : A 0 A* —>■ K is the cost function associated with the optimal control 
problem. 

Therefore, the optimal control problem consists on finding an admissible curve 
f{f) = y°'{t)) solution of the controlled Euler-Lagrange equations, the bound¬ 

ary conditions and minimizing the cost functional for C : A(B A* —>■ K. This optimal 
control problem can be equivalently solved as a second-order variational problem 
by defining the second-order Lagrangian L : A^‘^'> —>■ M as 


Liq\y‘^,r) = C{q\y^, 


dt 





(19) 


Here A^^^ denotes the set of admissible elements of the Lie algebroid A, a subset of 
A X TA, given by 

A^'^'> := {{b,Vb) G AxTA\ p{b) = TTA{vb)} 

where Tta '■ TA ^ TQ is the tangent map of the bundle projection. A^'^^ is con¬ 
sidered as the substitute of the second-order tangent bundle in classical mechanics 

[46] . In local coordinates, the set A^'^'> is characterized by the tuple (g*, ?/“, z“, u“) G 

A X TA such that ?/“ = z“. Therefore one can consider local coordinates (g*, y“, u“) 
on A^^K ^ 

The dynamics associated with the second-order Lagrangian L : A^^'^ —>■ K (and 
therefore the optimality conditions for the optimal control problem) is given by the 
second-order Euler-Lagrange equations on Lie algebroids (see for example [12] and 

[47] ) 


0 = 


d^ 

dt"^ 



■Clp{q)y^ 


d 

dt 



d dL 
Jt^ 


Clp{y)y^ 


dL 

dTf 



( 20 ) 


dq^ 

together with the admissibility condition —— = 

dt 
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Remark 6. Alternatively, one can define the Lagrangian L : — >■ K in terms of 

the Euler-Lagrange operator as 

L = Co 0 £C{L)) : ^ M, 

where £C{L) : A^^^ —>■ A* is the Euler-Lagrange operator which locally reads as 


£C{L) = 


/ d dL 

ydtdy^^ 



Here {e“} is the dual basis of {ca}, the basis of sections of A and : A^^^ —>• A is 
the canonical projection between A*-^^ and A given by the map A*-^^ ^ (g®, j/“, u“) i—)■ 
(<?*, 2 /“)eA o 


3.3.2. Optimal control problems of underactuated mechanical systems on Lie alge- 
broids. Now, suppose that our mechanical control system is underactuated, that is, 
the number of control inputs is less than the dimension of the configuration space. 
The class of underactuated mechanical systems are abundant in real life for differ¬ 
ent reasons; for instance, as a result of design choices motivated by the search of 
less cost engineering devices or as a result of a failure regime in fully actuated me¬ 
chanical systems. Underactuated systems include spacecrafts, underwater vehicles, 
mobile robots, helicopters, wheeled vehicles and underactuated manipulators. We 
will see that the corresponding optimality conditions are given by the solutions of 
second-order constrained Euler-Lagrange equations (see [16]). 

Given a Lagrangian function L : A —>■ M and control external forces, the con¬ 
trolled equations for an underactuated system defined on a Lie algebroid are 


^ ( dL\ 
dt \dy°-) 

dt \dy^ ) 


p: 


dL ^ g dL 

rg^i+ 




= {uF)a, 

= 0 , 


dt 


=Pay > 


with {a} = {A, a}. The optimal control problem consists on finding an admissible 
trajectory ^(t) = (g®(t), y^{t), y°'{t)) solution of the controlled Euler-Lagrange equa¬ 
tions given boundary conditions and minimizing a cost functional C : A 0 A* —> M. 

This optimal control problem can be solved as a constrained second-order varia¬ 
tional problem on Lie algebroids where the second-order Lagrangian L : A*^^^ —K 
is given by 


L{q\y^ 


)=C(g®,,^- 




( 21 ) 


and where the dynamics is restricted by the second order constraints 




f)T 


The optimality conditions for the optimal control problem are determined by 
the second-order constrained Euler-Lagrange equations given by considering the 
extended Lagrangian L = L-\- : A^^^ x K® —M where = (Ai,..., A^) S M® 

are the Lagrange multipliers. These equations are given by (see [12] for more details) 
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0 = 


together with the admissibility condition g® = 


3.3.3. Optimal control problems on Lie groupoids. Now we describe the discrete 
optimal control problem on a Lie groupoid G. Let : G —>■ K be a discrete La- 
grangian, an approximation of the action corresponding to a continuous Lagrangian 
L : A defined on a Lie algebroid A, that is, 

j■{k+l)h 

Ld(fffc) ^ / L{^{t))dt 

Jkh 

where h > 0 is the time step with T = Nh and ^ is an admissible curve on A. 

The discrete controlled Euler-Lagrange equations are 

- {rg^+, o i)*dhd{gk+i) = u/c € (22) 

for fc = 1 ,..., iV — 1 where go and g^ are fixed on G. 

We define the subset Gjs A*G of G x A*G, 

Gp A*G := {{g,u) &GxA*G\ P{g) = ta^g{u)}. 

Given a discrete cost function G^ : Gp A*G —>■ K, the discrete optimal control 
problem is determined by extremizing the discrete cost functional 

7V-1 

Jdigo^gi, ■ ■ ■ ,gN) ■= ^ Gd{gk,Uk), (23) 

k=0 

for (go, 5 i, • ■ ■ , 5 m) € G^+\ satisfying equations ( 22 ) with {gk,gk+i) € G 2 , fc = 
0,..., TV — 1 and with go, gi, gN-ij gN and g = gogi... g^ & G fixed points on G. 
Defining the discrete second-order Lagrangian : G 2 —t K as 

Ld{gk, gk+i) ■■= Gd {gk, tg^dLd{gk) - O i)*dLd( 5 fc-Hi)) , (24) 

the discrete optimal control problem consists on finding a path (go, 5 i, • ■ ■, 5 iv) € 

qN+1 

minimizing the discrete action sum J7d for the discrete second-order La¬ 
grangian Ld : G 2 -t M where go,5i,5 m-i,5m and 5 = 5051 ...gN € G are fixed 
points in G. 

Discrete Hamilton’s principle (3.2) states that the paths minimizing subject 
fixed points 50 ,5i,5 m-i,5m S G satify the discrete second-order Euler-Lagrange 
equations for : G 2 —t K given by 

{DiLd{gk,gk+i) + L'2Ld(5fc-i,5fc)) 

+ (Tg ;,+1 oi)* {DiLd{gk+i,gk+ 2 ) + D 2 Ld{gk,gk+i)) ■ 

Therefore, as in the continuous problem, the optimality conditions for the dis¬ 
crete optimal control problem are determined by the discrete second-order Euler- 
Lagrange equations for Ld : G 2 ^ R. 

Alternatively, one can start with a continuous optimal control problem associ¬ 
ated to a Lagrangian L : A —J, R defined on a Lie algebroid A. The optimality 
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conditions for this optimal control problem are determined by a system of fourth 
order differential equations obtained from the second-order Euler-Lagrange equa¬ 
tions associated with the Lagrangian L : —>■ M determined by the cost function 

as in (21). Now, we take directly a discretization of the second-order Lagrangian L 
to derive : G 2 —K. 

Finally, we would like to point out that the underactuated case follows as in 
the continuous case by consider a discrete second-order constrained problem as 
in Subsection 3.2 and the optimality conditions are given by the solutions of the 
discrete second-order constrained Euler-Lagrange equations (18). We will illustrate 
this in Example 3.3.5. 

3.3.4. An illustrative example: Optimal control of a rigid body on SO(3). In this 
example, we show how the optimal control problem of a rigid body defined on the 
Lie group SO{3) can be studied using the previous constructions given before. This 
example is motivated by the attitude optimal control of spacecrafts (see [35], [36] 
and references therein). 

The Lie groupoid structure of SO{S) over the identity matrix Id is given by 

a{R) = Id, /3{R) = Id, e{Id) = Id, i{R) = R~^ and m{RG) = RG 

for R,G £ SO{3). The Lie algebroid associated with the Lie groupoid 5'0(3) is the 
Lie algebra so(3) over a single point, where the anchor map is zero and the bracket 
is the usual commutator of matrices and the set of admissible elements is identified 
with so(3) X so(3). Observe that, in this case, all the elements are composable, that 
is, S 0 { 3)2 = 50(3) X 50(3). 

The equations of motion of the controlled rigid body are 

ill = 5*10203 -|- Ui, O2 = 520i03 -|- U2, O3 = 5*30102 -|- M3, ( 25 ) 

where O = (Oi, O 2 , O 3 ) S and O = (Oi, O 2 , O 3 ) £ Ui are the control inputs 
(torques for the rigid body) for * = 1, 2, 3 and 



are constants determined by the moments of inertia of the rigid body /i, 52 ,/ 3 . 
Here, we are using the typical identification of the Lie algebra so (3) with by the 
hat map t 50 ( 3 ) (^303 [2] and [29] for example), where with some abuse of 

notation, we directly identify with so(3) by omitting the hat notation. 

Our fixed boundary conditions for the optimal control problem are (i?(0), H(0)) 
and (i?(T), H(5)), where R{t) £ 50(3) is the attitude of the rigid body subject to 
the reconstruction equation R = RQ and variations for the attitude are given by 
SR = Rr], with rj an arbitrary curve on so(3). Consider the cost functional 



From eqs. (25) we can work out mi, M 2 and M 3 in terms of H and H. Consequently, 
we can define the function 5 : so(3) x so(3) — >■ K by L{0,tl) = ■ m(H,H) 

where u = {ui,U 2 ,u^). Therefore, the Lagrangian function L : so(3) x so(3) —>• K is 

5(H, H) = - ~ 5 * 1 ^ 2 ^ 13 ^ -|- ^^2 — 5*2HiH3^ -|- ^^3 — P3HiH2^ ^ . (26) 
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From (26) the cost functional becomes into J = 


L(n, tl) dt, (see [16] for the 


solution of this second-order variational problem in the continuous setting). 

Angular velocities and angular accelerations can be approximated by discrete 
trajectories ^k+i /2 — ^{kh) and ^k,k+i — ^(kh) respectively for ^k^^k+i G so(3) 
where h > 0 is a fixed real number and T = kh with k = 0,..., N, where we are 
using the notation ^k+i /2 = ^i^k + ik+i) and = j^i^k+i - &)• Define the 

second-order Lagrangian C : so(3) x so(3) K by £{^k,^k+i) = L {^k+i/ 2 ,ik,k+i) ■ 

The optimal control problem, is given by minimizing the cost function associated 
with the discrete second-order Lagrangian : S'0(3) x 5'0(3) —?► K over discrete 
paths on S'0(3) where 


Ld{wk,Wk+i) = hC{cSiy ^(wfc),cay ^(?«fc_|_i)). (27) 


Here {wk,Wk+i) G SO{2>) x SO{Z), c&y[h^k) = Wk and cay : so(3) —)■ 5'0(3) denote 
the Cayley map for the Lie group 5'0(3) (see Appendix D). 

Therefore, the discrete Lagrangian is now given by 

T , . , 7 /cay“i(w;fc) + cay“i(wfc-^i) cay“i(w;fc-n) - cay“i(wfc)\ 

Uiwk.Wk^,) = hL - - -,-- j . 

The variational integrator for the optimal control problem is given by applying 
discrete Hamilton’s principle (3.2) in the discrete action sum determined by the 
discrete cost 

N-l 

Cd = ^ Ld{wk,Wk+i). (28) 

k=0 

Instead of working with the discrete sum (28) one can take 


N-l 

Cd= L {Ck+l/2^^k,k+l) 
k=0 


(29) 


in order to work in a vector space where variations of ^k = cay ^{wk) G so(3) are 
given by (see for example [34]) 


S^k = ^ (^Ad^^r]k+i -r]k + ^ad^^r^fc - ^adj, {Ad^^r]k+i) + ^^kVk^k 


^ ik{Adwkdk+l)^k 


(30) 


3.3.5. Example: optimal control of a heavy top with two internal rotors. The fol¬ 
lowing example illustrates the study of underactuated mechanical control systems 
on Lie algebroids and the construction of variational integrators for such systems. 
It is the optimal control problem of the upright spinning of the heavy top (see [II] 
and reference therein). 


Consider the top with two rotors so that each rotor’s rotation axis is parallel to 
the first and the second principal axes of the top as in Figure I. Let /i, / 2 , /a be the 
moments of inertia of the top in the body fixed frame. Let Ji, J 2 be the moments 
of inertia of the rotors around their rotation axes and Jji, Jj 2 , Jjs be the moments 
of inertia of the j-th rotor, with j = 1,2, around the first, the second and the 
third principal axes, respectively. Also we define the quantities /i = /i -I- Jn -I- J 21 , 
I 2 = I 2 + J 12 + T 22 , 13=13 + Ji3 + <^23) 7 i = -I- Ji and 72 = 72 + J 2 - 
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n 


Center of Mass 


Fixed point ^ 


Figure 1. Heavy top with two rotors, each consisting of two 
rigidly coupled disks. 


Let M be the total mass of the system, g the magnitude of the gravitational 
acceleration and h the distance from the origin to the center of mass of the system. 

The system is modeled on the transformation Lie algebroid H x so(3) x 

T(S^ X S^) over the manifold Q = x (S^ x§^), where the anchor map p : 
X so(3) X r(S^ X S^) —>■ x x S^) is locally given by 



Here H = (Hi, H 2 , H 3 ) G so(3) ~ is the angular velocity of the top in the body 
fixed frame, L = (ri,r 2 ,r 3 ) G represents the unit vector with the direction 

opposite to the gravity as seen from the body and 0 = ( 0 i, 6 * 2 ) is the rotation angle 
of rotors around their axes. 

If we denote by (z = 1, 2, 3), the standard basis of matrices of so(3), given by 



0 0 0 

0 0-1 
0 1 0 



0-10 
1 0 0 

0 0 0 


then the basis of sections of A is given by the elements 



with 1 = 1,2, 3. 
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The Lie bracket of sections of A is determined by 

= X^^, 

_ j^[Ei,E3] _ _^E2^ 

with r, s = 1 , 2 and , X^'j = 0 , for s = 1 , 2 and z = 1 , 2 ,3. 

The reduced Lagrangian ? : x so(3) x r(S^ x S^) —>■ M is given by 


i(r,o,0) 


/ \ 

T 

0 

0 

0 


/ \ 

il2 


0 72 0 0 J2 


fl2 

ils 


0 

0 

0 

0 


Q3 

01 


0 

0 

0 


01 

[ 02 ) 


0 

0 

0 


\ 02 / 


— MghT^. 


The Euler-Lagrange equations for I are given by 



/ dl 

dt 



( dl 

dt ' 



= bd. ^ ^ + MghT X 63, 

oil 

= 0 , i = l, 2 ; 


together with the admissibility condition f = T x fl. 

Next, we add controls in our picture. Each rotor can be controlled is such way 
the controlled Euler-Lagrange equations are now 


d 

/ dl 

dt 


h 

( dl 

dt ' 



dl 

= ^ X fi + MghT X 63, 

—ZZZ, Z — I5 2 , 

r =r X 


where 63 = (0,0,1). That is, 

7i^i T ~ 72 ^ 2^3 + ^ 313^2 =MghT2, 

72^2 + <^2^2 + 71 ^ 1^3 ~ >.^i^if^3 = ~ AlghTi, 

^3^3 ~ 71 ^ 1^2 ~ Jl0lO2 + 72^12^^! + J2^2fll =0, 

Ji(f 2 i + 6 i) =ui, 

J2{^2 + ^ 2 ) =U2, 

together with the admissibility conditions 

Li = r2n3 — r3n2, r2 = r3r2i — ri03, r3 = rin2 — r20i 
where T = (Li, r 2 , r 3 ) G 

The optimal control problem consists on finding an admissible curve 7 (t) = 
(r(t),fl(t),(l(t),Ui(t)) of the state variables and control inputs, satisfying the con¬ 
trolled equations given above, the boundary conditions and minimizing the cost 
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functional 

1 r'^ 

= - y (ui + U 2 ) dt. 

This optimal control problem is equivalent to solve the second-order variational 
problem determined by 

min f L dt, 

and subjected to the second-order constraints : IR^xIR^x2so(3)xT*^^^(S^xS^) —)■ 
M, A = 

0, r, n, 0 , t, 0 ) = 7 iiii -t Jl01 - - MghT 2 = 0, 

0, r, O, 0, r, = 72^2 + •/202 + 71 ^ 1^3 ~ + MghT 1=0, 

0 , r, o, 0 , r, 0 'j =i^ri 3 — 7 ini 02 — Ji 0 iii 2 + 72 ^ 2 ^^! + J 202 ^i = 0 , 

together with the admissibility condition F — T x = 0 and where L : so (3) x 
so(3) X X S^) —)• K is given by 

L(n,0,ii,0,0)=c(n,0,n,0,f^(^^"j^ = :^{n, + 0,f+ :^{ii2 + 02f 

where C = ^(ui -I- m|). 

Therefore, the optimality conditionsare determined by the constrained second- 
order Euler-Lagrange equations given by 

0 + ^1 ^) + Ai — A2i23 -|- A2fl3 — A3ii2) 

0 =J2{’ii2 + ^2 ^) + A 2 — A3J7i — A 2 AI 1 , 

0 =7iJ7i -|- — 72^12^3 -f $ 23/3122 ~ MghT2, 

0 =72$22 -f J202 + 7 iA^i^3 ~ >2101123 -f MghT 1 , 

0 — 7i12i122 — Ti0i122 -I- 72 I 22 I 21 -I- J202^ii 

0 =A{X^B A 3 ) -f A 3 xi -|- Ji ('$2i -f 01 ^ — 12ii? -|- 0 i_B) 

T A27i(^3 T 12i123 -|- I 22 I 23 -|- 123 ) -|- 7i(A 123 — Ai — Ai-B), 

0 =AiC -f (A 3 -f Ai?)(72l2i — 7 i12i — J\0i) — A 272 + AiC* 

~ A3(72$2i — 7i$2i — Ji0i) — J|(122 -I- 02 ^ -I- -B) -I- B{\iC — 72 A 2 ) 

0 =(A 2 + -BA2)(7 i12i — Ji0i) — B\i{'^2^2 + 23 I 22 ) — BI^{\iTl 2 + A 3 ) 

+ A 2 ( 7 i^i ~ < 2 i 0 i) — 72(AiAl2 + Ai 122) — /3(Ail22 — A3 — 2Ai$22 — Ail22 

0 =f - r X 12 

where A = 72 II 2 - 7 i $22 -f J 202 , i? = 12i -I- 122 -I- I 23 and C = 113/3 - 72 ll 3 - 

As before, we use the Cayley transformation on SO{3) to describe the discretiza¬ 
tion of the optimal control problem for the heavy top with internal rotors. We 
redefine the Lagrangian L and the constraints as C : 2so(3) x x S^) —>■ M 
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and : TR^ x 2so(3) x r(2)(§i x S^) ^ M by 

^{^k,&,^k+i,&, 0 ) := L (4+1/2, ^,C/c,/c+i, 


and 


$^(r,4, 0, f,4+1,0, 9 ) := (r,4+1/2, 9, f,4.fc+i,0), 

where 4 , 4 +i ^ so(3) and h > 0 is a fixed real number with T = kh, k = 0,..., N. 

The discrete second-order Lagrangian Ld : 2so(3) x 3(S^ x S^) —)• M associated 
with C is given by 


-bd(4, 4 , 4+1, ^fc+l, ^fe+ 2 ) ■“ ( 4 + 1 / 2 , ^fc+ 2 / 3 , ^fc.fc+l, ^fe,/c+2, ^fc+2//i2 


for z = 1,2 and the discrete constraints : 2so(3) x 3(S^ x S^) —)■ M associated 
with by 


(Ffc, 4, ei, 4 + 1 , Tfe+i, 0 ^ 1 , 0 ^+ 2 ) := 

(Pfe, 4 + 1 / 2 , (^k+2/3’ rfc,fc+l, 4,fe+l, 4,fc+2, ^k+2lh^ 


for z = 1, 2 , where 


0 


fc+2/3 


+ ^fe+l + ^fc+2 


, 0\ 


k,k+2 


ni ni 

^fc+2 


91 


2 h ’ '=+2/'^^ 


ni nni i ni 

^fc+2 


and 0-4 = cay ^{ujk) G so(3) for ujk € SO{3). Here 

0 -(4)fe (4)fc 

4=1 (4)fc 0 -(4)fc 1 gso(3). 

-(4)fc (4)fe 0 


The geometric integrator is given by extremizing the discrete cost function de¬ 
fined by 


N-l 


Cd=J2 [i44,0L4+i,0l+i,0U2) + (A^f'f>d(rfc,4, ^fc,rfc+i,4+i, ^fe+i, ^fc+ 2 )] 


k=0 


where A+ G A = 1,...,3 are Lagrange multipliers. That is, it is given by 
the solutions of the discrete second-order constrained Euler-Lagrange equations 
associated to the discrete extended Lagrangian Ld = Ld + ^A^d ■ 2oo(3) x 3(S^ x 
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T '^1 /ir I I \ AC/c+i)2 - (^fc)2 (^k+2 

Ld = Y + Ou+2/h^) + Y - h -+- W~ 


q2 ^ 2 


(j)l ~ I 7. /^^fc + 2 “ ^^fe + l + 


■ 72 

■ h 


h J ' 

^ ((C/c)2 + (^fc+l)2)((C/c)3 + (Cfc+l)3) 
((?fe+l)3 — (C/c)3)((?fe+l)2 + (C/c)2) 


- MghYk 


2 h 


»3 I KiiilY^ )+*'"*+" 


n2 0/32 I n2 


+ 7i 
— Ji 


h 7 ' V 

^ ((C/c)l + (^fc+l)l)((gfc)3 + (Cfc+l)3) 

i^k+2 ~ ^k)ii^k+l)3 + (Cfc)3) 


2 h 


- MghYk, 


s3_f 7(Cfe+l)3 - (^fc)3^ 7((Cfc+l)l + (?fe)l)((?fe+l)2 + (C/c)2) 


l&d =I3 


V h j 

^ ((C/c+l)2 + (^fc)2)((C/c+l)l + (^fc)l) 
(^fe+2 ~ ^D((?fc)l + (Cfc+l)l) 


■72 

■ -h 


Ah 


— Jl 


i^k+2 (^k)ii^k)2 + (Cfc+l)2) 


together with Ffe+i = 


Ah 

cay~^{ujk) + cay~i(a;fc+i) 
2 h 


A. Lagrangian submanifolds generating discrete dynamics. In this section 
we study how a Lagrangian submanifold of a particular cotangent groupoid can be 
used to give a more geometric and intrinsic point of view of discrete second-order 
dynamics. Moreover, we will study the preservation properties of the derived dis¬ 
crete implicit dynamics. We also study discrete second-order constrained systems. 
Particularly, from this geometrical framework, we will analyze some geometric prop¬ 
erties of the associated discrete flow. Finally, we will study the theory of reduction 
under symmetries. 

4.F The prolongation of a Lie groupoid over a fibration. Given a Lie 
groupoid G =1 (5 with structural maps a:G—>-(5, /3:G—J-Q, eiQ—>G, 
i : G —>■ G, m : G 2 —>■ G, and a fibration tt : P —)• Q, we consider the set 


V^G = Xa Gpx^P = {{p,g,p') € P xG X P/tt{p) = a{g),l3{g) = tt{p')}. 
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V^G has a Lie groupoid structure over P, where the structural maps are given by 
-.r^G—> P, ip,g,p')< — >p; 

irp : {V^G )2 — {{p, g,p'), {p', h,p”)) i—s> {p,gh,p”); 

■. P — >V^G, pi —> {p,e{a{p)),p)] 
r:p-G^P"G, {p^g^p')^{p',g-\p). 

V^G is called prolongation of G over tt : P ^ Q (See [40],[41] and [56]). 

Next, we consider the prolongation V°‘G of the Lie groupoid G over its source 
map a : G ^ Q, that is, one can consider the subset of 3G := G x G x G, 

P“G = G„ x„ G^ x„ G = {(p, h, r) e 3G/a(5) = a(h), /3(h) = a(r)}. 

P“G is a Lie groupoid over G. Moreover, G2 Q P“G where the inclusion is given 

by 

^G, : G2 P“G 

(g,h) I—s> (g,g,h). 

Now, we construct the Lie algebroid associated with P“G. This will be identified 
with the prolongation P°‘(AG) of AG over a : G Q, where AG is the Lie algebroid 
associated with G with bundle projection tag '■ ^G —> Q. 

Definition 4.1. The Lie algebroid associated with a prolongation of a Lie groupoid 
G over a is given by, 

Ag (P“G) = {(a,(„(g));rg) e A„(g)G x T,G | (T,a)(Tg) = (r.(a(,))/3)(a,(„(g)))}. 

A (P“G) is a Lie algebroid over V’^G with bundle projection denoted by ta{V‘>g) '■ 
A(P“G) ^ P“G. 

Remark 7. The prolongation of a Lie algebroid over a fibration has a Lie algebroid 
structure (See Appendix B). If we consider the linear isomorphism 

(vl/“)g : A, (P“G) ^ P“(AG) C A,(g)G X T,G 

(Og, Og((j(g)) , Ij)) I > (®£(a(g))l bg) , Vp G G 

then the mapping (4'“)g, g € G induces an isomorphism ilt" : A (P“G) —>■ V°‘[AG) 
between the Lie algebroids A (P“G) and V^^AG) (see [28] and [41] for more details 
about the construction of this linear isomorphism). o 

From (4.1) a section Z oi A (P“G) can be expressed as 

Z{g) = {X{a{g)),Y{g)) 

where g G G, X G r(TAG) and T is a vector field on G such that TI3{X) = Ta{Y). 

The corresponding left-invariant and right-invariant vector fields associated with 
the section Z G ^{tacp<^g)) are 

Z(g,/i,r)=(0,^(h),y(r)), (31) 

t{g,h,r)={-Y{g),l{h),0), (32) 

with {g, h, r) G V^G. 

Given a basis of sections of AG one can obtain a basis of sections of A (P“G), 
denoted by {Zi, Z 2 }, with 

Zi = i-X,^) and Z2 = (0,^) 


(33) 


24 


LEONARDO COLOMBO AND DAVID MARTIN DE DIEGO 


where X G A G a(G) and a G X(G). Here, we are using the notation 

1^(G) (resp., X (G)) for the set of right-invariant (resp., left-invariant) vector fields 
on G. 

The next result is a direct application of the construction given before and it 
will be useful when we derive the discrete second-order dynamics for a second-order 
discrete Lagrangian. 

Lemma 4.2. Let {Zi,Z 2 } be a basis of sections of AlV^G), where Z\ = (—X, 
and Z 2 = (0, X G r(r^G)- For {g,h,r) G 7^“G the associated left and right 
invariant vector fields for Zi and Z 2 are given by 

^i{ 9 ,h,r) = (-^(g),-]?(/i),0) ^ 2 ( 9 , h,r) = (-^(g),0,0) 

%{ 9 ,Kr) = (0,-^(/i),^(r)) %{ 9 ,h,r) = (0,0,^(r)). 

4.2. Generating Lagrangian submanifolds and dynamics on Lie groupoids. 

Let G =} Q be a Lie groupoid with source and target map a, ft : G —)■ Q respec¬ 
tively, and we consider the prolongation of G over its source map, 'P°‘G. We denote 
by a“,/3“ : 7^“G —>■ G the source and target maps of this Lie groupoid. Let 
ta*{V‘^G) ■ A*{fP°‘G) —>■ G be the dual of the vector bundle associated with the 
Lie algebroid ■ A{V°‘G) -A G. The Lie groupoid (cotangent groupoid) 

T*{V°‘G) =t A*{V°‘G) is a symplectic groupoid (see example 3 in section 2.1.1). 

In what follows, we show how the discrete dynamics associated with a discrete 
second-order Lagrangian : G2 —>■ K is generated by a Lagrangian submanifold of 
the cotangent groupoid T*{V°‘G). 

Remember that given a manifold Q and a function S' : Q —>■ K, the submanifold 
Im dS C T*Q is Lagrangian. There is a more general construction given to Sniatycki 
and Tulczyjew [57] (see also [58] and [59]) which we will use to generate the discrete 
dynamics. 

Theorem 4.3 (Sniatycki and Tulczyjew [57]). Let Q be a smooth manifold, N C Q 
a submanifold, and S: iV —>■ K. Then 

Ss = {^ G T*Q 1 TTQ{fi) G N and {g,,v) = {dS,v) 

for all V G TN C TQ such that tq{v) = 7rQ(/i)} 

is a Lagrangian submanifold ofT*Q. Here ttq : T* Q ^ Q and tq : TQ ^ Q denote 
the cotangent and tangent bundle projections, respectively. 

Turning back to the groupoid formulation, immediately from Theorem (4.3), the 
discrete second-order Lagrangian : G2 —>■ K generates a Lagrangian submanifold 
C T*{fP°‘G) of the symplectic Lie groupoid {T*{fP°‘G),uj'pcQ) where w-pcQ 
denotes the canonical symplectic 2-form on T*(7^“G). Denoting by ^^2 ■ G 2 ^ 7^“G 
the inclusion defined by *02(51152) = (91,9i, 92), we have 

Sl. = {5 e T*{r^G) 1 = dLd} C T*(P“G) 

is a Lagrangian submanifold of {T*ifP°‘G),oj'P‘^G)- 

The relationship among these spaces is summarized in the following diagram 
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c-^ A*{V°‘G) 

p 

T ir. » a“ '' 

K G 2 ^^ ? G 

/3“ 

where from now on we will denote a and /3 the source and target maps, respectively, 
of the Lie groupoid T*{V°'G) =1 A*(V°'G). 

Given an element ^ S ^ (7^“G) with {g, h, r) S 7^“G the source and target 

maps of T*{V°'G) are defined such that, for all section Z G ^{tacp^G))^ 

(5(/r), Z{a{g))) ={g., 't {g, h, r)) (34) 

0{fA,Z{P{g))) ={g,^ig,h,r)), (35) 

where Z and ^ are the corresponding left and right invariant vector helds associ¬ 
ated with the section Z of A{V°‘G) according to (31) and (32). 

Denoting by 

^(gk,gk+i) = (Msfe)Msfc)i*gfc+i) € ^icsGfc.gfc+il^'^ 

with {gk,gk+i) € G2, the Lagrangian submanifold gives rise to the discrete 
second-order dynamics as we describe in the following. 

Definition 4.4. A sequence 7(gi,g2 )> • • • iTIsn-i.sn) ^ T*{V°‘G) satisfy the second- 
order dynamics on if 7(gi,s2)> ■ • ■, 7 (sjv-i, 3 n) ^ and 

^ilGk.gk+i)) = hl{gk-i,gk)) for fc = 2,..., iV - 1. 

That is, 7(51,92)1 ■ ■ ■:7(giv-i.gjv) composable sequences on T*{V°‘G). 

Theorem 4.5. Let G ^ Q be a Lie groupoid and L^ : G 2 —>■ K 6 e a discrete 
second-order Lagrangian. Consider the Lagrangian submanifold of the cotan¬ 
gent groupoid T*{V°‘G) generated by Ld- A segttence 7 ( 91 , 92 ), ,7(gjv-i,gw) satisfies 

the discrete second-order dynamics on if and only */ 7 ( 91 , 92 ), ■ • •, 7 (gjv-i,giv) 
isfy 

i9k)j P-gk) “b O^idk)^ Pgk+i) idk+l) j Pgk+i) ~ 0^ idk-i-l) ^ h-gk+i) 

Lgk “b Pgk ~^l^digkj gk-kl) 

Pgk+i =D 2 Ld{gk,gk+i) 

for k = 1,N — 1 and for any section Z G L{ta(V'^g)) according to (31) and 
(32). 

Proof. Consider the sequence 7(91,92), • ■ ■ , 7 (gjv-i,gN) fo T*{V°‘G). Applying the 
definition of a (34) and /3 (35) to the relation <5(7(9,,,91,^1)) “ Pil{gk-i,gk)) fo^ 
k = 2, ...,A^ — 1, we have that for any section Z G L{TA(j>o,Qf) the sequence 
7(91,92)^ ■ • ■ ,7(gjv-i,gN) belongs to the Lagrangian submanifold if 

( ^ i.gki gk^ 5fc+l), 7(gfc,gfc+i)) ("^(Sfc+l, gk+h <7fc+2), 7(gfei-i,gfc+2)) 

Lgk “b Pgk —DiLd^gki gk+i) (3f) 

Pgk+1 =D2Ld{gk,gk+i), (38) 
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for /c = 1,..., — 1. Using (31) and (32) the equations given above are equivalent 

to 

0^i9k)i P'Qk) "b (^(fffe + l)) Mgic+i) {9k+l)i P-gk+i) ~ )) A^gt+i) (39) 

Nk + 9-gk =Di Ld i9k ,9k+i) (40) 

Asfe+i =D2Ld{gk,9k+i) (41) 

for fc = 1,..., -/V — 1 as we claimed. □ 

Remark 8. We have seen how the dynamics is only defined implicitly by a relation 
in T*{V°‘G) rather that as an explicit discrete flow map. Therefore, the sequence 
7(9 i,S 2)’■ • ■ ’7(giv-i.gN) ^ T*{V°‘G) satisfy the discrete second-order dynamics on 
Ez, if and only if each pair of successive elements satisfies the relation 

(7(gfc_i,gfc),7(gfc.gfc+i)) ^ {T* {V°‘G))2 H (El^ X E^J. 

o 


Next, we show that the discrete dynamics described implicitly in Theorem (4.5) 
is equivalent to the discrete second-order Euler-Lagrange equations (15) given by 
the variational point of view. 


Theorem 4.6. Let : G 2 —> K 6e a discrete second order Lagrangian. For 
every section Z ofV^TAijicQ^) according to (31) and (32) the discrete second order 
Euler-Lagrange equations associated with L^ are 

^ ~^gk {DiLd{gk,gk-i-i) + L)2Ld{gk-i,gk)) 

+ {’’^gk+i ° 0 * (DiLdigk-kii gk+ 2 ) + D2Ldigk,gk-i-l)) ) 
for k = 2,...,N -2. 


Proof. Let Z be a section of A{V°'G) and consider the basis {^ 1 ,^ 2 } of section 
of A{V°‘G) as in (33). Using Lemma (4.2) in (36) we get the relation 

(^2 (.gk j gk^ l?fc-t-l); ^{gk ,gfe+i)) (^2 (l/fc-t -1 : ; I//C+2) 1 7 (gfc+i ,gfe+2)) 


if and only if 

0^{9k+l) j9gk+i) — (^(fffeH-l)) Mgfc + i) (42) 

for fc = 1,..., IV — 1. Now, using the relations 


h-gk+i 4 " Lgk+i —Di^d{gk-i-ij gk+2), 
Ugk+i —Gt2Ld{gk-, gk+i)j 


(43) 

(44) 


we have, 

(0 (gk+i),D2Ld{gk,gk-ki)) = - {0 {gk-ki),DiLd{gk-i-i,gk-i-2) - Agt+i)- (45) 
Similarly, by Lemma (4.2) in (36), we get the relation 

( ^ l(l7fe) 5fcj 5fc-l-l)) 7(gfc,gfc + i)) ("^1 (Hfc-l-lI 5fc+l) 5 ^+ 2 )) 7(gfc+i,gfc+ 2 )) 

if and only if 

(0 {9k), idgk) ~ (0 {9k+l) J Lgk+l ) — (^(5fc-i-i) ! Lgk+1 ) + (^(5fe-6l),Agfc+i)- (46) 

for fc = 1,..., iV — 1. Using (43), (44) and (45) equations (46) are equivalent to 

0 ={0{9k-i-i), L)iLd{gk-i-i, 9k+2) + D2Ld{9k, 9k+i)) 

— {gk), DiLd{gk, 9k+i) + D2Ld{gk-i,gk)) 
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i.e., 

~^*gk {DiLdigk-, Qk+i) + D2Ld{gk-i,gk)) 

+ (^sfc+i ° *)* {DiLd{gk+i, gk+ 2 ) + D 2 Ld{gk, gk+i)) 
for k = 2,..., N — 2, after a shifting of the indexes, as we claimed. □ 

Example 1. Let G be a Lie group and let : G 2 = G x G —)■ M be a discrete 
second-order Lagrangian. The prolongation of G over its source map, V^G, is a Lie 
groupoid over G and it can be identified with three copies of G, that is, 'P“G ~ 3G. 
We construct a Lagrangian submanifold of the cotangent groupoid T*{V°‘G) as 

Sl, = {/X e T*(iP“G) I = dLd} C r*(iP“G) ( 47 ) 

where i-p‘^G ■ G 2 —>■ 'P“G is the inclusion given by iv’^cigih) = {g,g,h) with 

{g,h)&G2. 

Observe that *X{gk) = Tig^{X), J^igk) = -T{rg^oi){X) = Trg^(X). Therefore, 
we have that a sequence where 7 (g,, 3,_^7 e T*^^^g„g,^,)i'P°‘G), 

with {gk, gk+i) € G 2 for fc = 1,..., — 1, satishes the discrete second-order dy¬ 

namics on if 

0 =^gkddiLd{gk,gk+i) + tg^D 2 Ld{gk-i,gk) - r*g^^^D 2 Ld{gk-, gk+i) 

- ‘r*g^^+^DiLd{gk+i,gk+2), 

that is, {gk, gk+i) € G 2 for A: = 1,..., TV — 1 is a solution of the discrete second-order 
Euler-Poincare equations for the discrete second-order Lagrangian : G 2 —>■ ffi. 


Example 2. Let Q be a differentiable manifold, consider the pair gronpoid QxQ ^ 
Q, where the source and target maps are given by the projections onto the fist and 
second factor, respectively. The set of admissible elements is given by 

{Q X < 3)2 = {(( 9 o,<?i),( 9 i, 92 )) & {Q X Q) X {Q X Q) \ qi = qi} ~ 3Q. 

The prolongation Lie groupoid is a Lie groupoid over Q x Q given by 


'P“(Q X g) = {((go,<?i), ( 92 , 93 ), ( 94 , 95 )) e 3(g xQ)\qi=q 2 and qs = < 74 } 4g, 

where the inclusion of (g x g )2 into V°‘{Q x g) is given by 

* 3 Q : (g X g )2 77 3g iP“(gxg)~4g 

( 90 , 91 , 92 ) I—^ ( 90 , 91 , 91 , 92 )- 

Given Ld '■ {Q x Q )2 —)■ ffi, a discrete second-order Lagrangian, we construct the 
Lagrangian submanifold of the cotangent groupoid (T*(7^“(gxg)),w-po(QxQ)) 
over A*{V°'{Q x Q)), where w-po^g^Q) denotes the canonical symplectic 2-form on 
T*(7^“(gxg)), by 


Sp, = {/I G T*(7^“(g X g)) I z*q 7 i = dLd}. 

Here n = fiodqo + fiidqi + gidqi + /i 2 d 92 - Therefore, /i if it satisfies 

^Ld , X _ dLd , X dLd . x^ 

Mo = ^—(90,91,92), Ml + Ml = ^—(90,91,92), M2 = ^—(90,91,92)- 
dqo dqi dq 2 

Using the source and target map given by 

d : T *{ V ^{ QxQ ))^ T *{ QxQ ) /ii, 7x2) ^ (-Mo,-Mi) 

/3 : T*{V^{QxQ))^T*{QxQ) (7Zo,Mi,Mi,M2) ^ (Mi,M2); 
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we have that the second-order discrete dynamics on holds if and only if 

D 2 Ldiqi,q 2 , qs) + DiLd{q2,q3, 94 ) + D^Ldiqo, qi, 92 ) = 0. 

4.3. Regularity conditions and Poisson structure. We have seen how the 
dynamics is implicitly defined by a relation on T*(P“G) rather than an explicitly 
defined map and pointed out that 7 ( 91 , 92)5 ■ • ■ )7(gN-i.gN) ^ T*{V°‘G) satisfies the 
discrete second-order dynamics if and only if for each pair of successive elements in 
r*(P“G) they satisfy 

(7(9.,9.+d5 7(9.+i.s.+.))e(r“'(P“G))2n(E,., xEiJ, A: = l,...,fV-2. (48) 

Weinstein [62] raised the question of how regularity results for the pair groupoid 
Q X Q might be generalized to arbitrary Lie groupoids G =4 Q, and this question 
was answered by Marrero et al. (see, Theorem 4.13 in [41]). Here, we study an 
extension of this problem to discrete second-order systems following the work of 
Marrero et al. [43] for first order systems. The problem consists on finding under 
which conditions the relation (48) is the graph of an explicit flow 

^(.9k — 1^9k'} ^ ^ ^(.9k T9k3-l') 

(at least locally) and what properties have such map. 

Consider the source map of the cotangent groupoid T*(P“G) restricted to the 
Lagrangian submanifold E^^^, that is, a ^Ld A*{V°‘G). If this map is a 
local diffeomorphism, then the Lagrangian flow is locally given by 

^Ld = (a ^Ld —^ ^Ld- 

Theorem 4.7 (Marrero, Martin de Diego and Stern [43]). Let T a symplectic 
groupoid over a manifold P with source and target maps a and /3, respectively. Let 
T, be a Lagrangian submanifold E C T. Then the restricted source map a |s: E —>■ P 
is a local diffeomorphism if and only if the restricted target map /3 |s: E —>■ P is a 
local diffeomorphism. 

A direct consequence of Theorem (4.7) is that given the symplectic groupoid 
{T* (V°‘G),ui'pog), if '^Ld C T*{V°‘G) is the Lagrangian submanifold generated by 
the second-order discrete Lagrangian : G 2 —t K, 5 '■'^Ld ^ A*{V°'G) is a lo¬ 
cal diffeomorphism if and only if /? : E^^^ —>■ A* (P“G) is a local diffeomorphism. 

The applications a [s^^ and jd [s^^ play the role of FLf and FPj, respectively, in 
discrete mechanics (see Appendix B and [13] for the higher-order case), that is, 

FL+=/3|e^_^ and FP-= <5 

and therefore, from Theorem (4.7), FP)]' is a local diffeomorphism if and only if 
^Lf is a local diffeomorphism. 

Now, by Definition (4.4), given a sequence 7 ( 91 . 92)5 5 7(g«-i,s«) ^ P*(P“G), 

it satisfies the discrete second-order dynamics on E^^j if and only if 7 ( 91 , 92 ); ; 

7(gjv-i,9iv) S '^Ld S-nd 

h{9k,9k+l)) = (7(g.+i.9.+2))5 k = l,...,N-2. 

Definition 4.8. Let P^j : G 2 —t K be a discrete second-order Lagrangian. It is said 
to be regular if a [s^^ is a local diffeomorphism, and hyperregular if it is a global 
diffeomorphism. 
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Definition 4.9. If the Lagrangian : G 2 —>■ K is hyperregular, the map defined as 

= P o(d : A*{V°‘G) —?> A*{V°‘G) is called Hamiltonian evolution 

operator. 

The next theorem shows the relation between the Hamiltonian evolution operator 
and the preservation of the Poisson structure on A*{V°'G). 

Theorem 4.10. Assume that a ® global dijfeomorphisms. Then, the discrete 

Hamiltonian evolution operator : A*(V°‘G) —>■ A*(V°‘G) preserves the Poisson 
structure on A*{V°‘G). 

Proof. If a (or, equivalently, (3 is a global diffeomorphism, the Hamilton¬ 
ian evolution operatior P^^ is a global automorphism on A*{V°‘G). 

Consider the application 

(a,/3) : T*{V°‘G) A*{r°‘G) x A*{r°‘G) 

T '—^ iHT),PiT)) 

where A*{V°'G) denotes A*{V°‘G) endowed with the linear Poisson structure changed 
of sign. 

The submanifold C T*{V°‘G) is the graph of f in A*{fP°^G) x A*{V°‘G). 
Since {a,/3) is a Poisson map and is a Lagrangian submanifold, the image 
of by {a, (3) is a coisotropic submanifold of A*{V°‘G) x A*{V°‘G). Thus, by 
corollary (2.2.3) in [61], Pl^^ is a (local) Poisson automorphism on A*{V°‘G) and 
therefore P^^ : A*{V°‘G) — >■ A*(V°‘G) preserves the linear Poisson structure on 
A*{P°‘G) as we claimed. □ 


4.4. Morphism, reduction and Noether symmetries. In this subsection we 
study the reduction of discrete second order Lagrangian systems and Noether sym¬ 
metries. Consider two Lie groupoids G =4 Q and G' Q' with structural maps 
denoted by a, 13, i, e, m and a', j3', i!, e', m' respectively. 

Definition 4.11. A smooth map % : G — >■ G' is a morphism of Lie groupoids if, 
for every composable pair {g,h) S G 2 , it satisfies (x(5),x(^)) G G 2 and xid^) = 
x{9)xW where G 2 denotes the set of composable pairs on G' Q' 

A morphism of Lie groupoids x : G —>■ G' induces a smooth map Xo ■ Q ^ Q' 
such that 

a' o X = Xo o a, o X = xo o /?, and x o e = e' o xo, 
that is, the following diagram is commutative. 


G 


a.,13 


^ Q 


X 

G' 


oc' ,0' 


XO 

I 




A morphism of Lie groupoids x induces a morphism Ax : AG —?► AG' of their 
corresponding Lie algebroids and 

%dv){x{9)) = TgxCv ( 5 )), 

Axiw]ix{g)) = Tgxi^ig)), 


(49) 

(50) 
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for all g S G, T € yl^(g)G and w S Moreover, for X € r(AG) and X' € 

r^AG') we have Ax o X = X' o xo and only IfTyo^ = X' ox (resp., Ty ojl = 
X'ox), (see [40] for more details). That is, X and X' are “Ay-related” if and only if 
their corresponding left-invariant (reps., right invariant) vector fields are y-related. 

Now consider the prolongations of G and G' by the source map a and a', respec¬ 
tively, denoted by 7^“G and V°‘ G'. 

Definition 4.12. Let y : 'P“G —>■ P" G' be a morphism of Lie groupoids and 
X G V^G. Two covectors ^ G T*{V°‘G) and G G') are said to be 

y*-related if (^,C) = {fJ-\TxiO) all ^ G Tj;{V°‘G). Also, if z G A*(7^“G) and 
^ ^ A*y-related if (z,5 = {z'^AxiO) for all f G A^(V°‘G), 

where yo : 'P°‘G —)■ G' denotes the smooth map on the base induced by the 
morphism y and where Ay : A{V°‘G) —>■ A{V°‘ G') is the associated Lie algebroid 
morphism. 

The following theorem states the reduction of discrete second order Lagrangian 
systems on Lie groupoids. It follows Theorem 4.5 in [43] for discrete first order 
systems. 

Theorem 4.13. Consider two Lie groupoids G ^ Q and G' ^ Q'. Let : G 2 —>■ 
K and L); : G 2 —>■ R &e a discrete second order Lagrangians, and let y : 7^“G —>■ 
G' be a morphism of Lie groupoids satisfying G 2 = y~^(G 2 ) and Ld = L'^ox 1g2 ■ 
If p G T*{V°‘G) and p' G T*{V°‘ G') are x* —related, the following properties hold 

1. If p' G then p G 

2. The sources a{p) G A*{V°‘G) and a'{p') G A*{V°‘'G') are A*x-related. 

3. The targets P{p) G A*(fP°‘G) and P'{p') G A*{V°‘ G') are A*x-related. 

Proof. (1) Consider p G Tf{V°‘G) and v G Tx{'P°‘G) with x G 7^“G. Since p and 
p' are y*—related we have 

{p,v) = {p',Tx{v)) = {dL'd,Tx |g 2 [v)) 

= {x*{dL'd),v) = {d{L'dOx Igs),'!') = {dLd,v) 
and therefore p G 

(2) Consider v G Aa(x)i'P°'G). Then {a{p),v} = {p,lt) = {p',Tx{1^)), because 
are y*—related. Using (50) we have 

{p',Tx{-^)) = {p\A^) = {A'{p'),Ax{v)). 

Therefore {a{p),v) = {a'{p'), Ax{v)), and thus a{p) and a'{p') are A*y- 
related. 

(3) Consider v G Tx{V°‘G) with x G 7^“G and observe that 

0{p),v) = (p,^) = {p',Tx(lI)) = {p',Ax(v)) = 0'{p'),Ax{v)). 

Thus, /3(/r) and P'{p') are A*y—related. 

□ 

Corollary 1. Let y : 7^“G —> 'P°‘ G' he a morphism of Lie groupoids. 

Ifl[g/ g')) • ■ • >7(g' g' ) ^ T*lfP°‘ G') sotisfy the discrete second-order dynamics 
for the discrete system determined by L'j^ : G 2 ^ ^ then any seguence x* — related 
7(91,92)’7(gN-i.sN) S T*{V°‘G) satisfy the discrete second-order dynamics for 
the discrete system determined 6j/ : G 2 —>■ K. 
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Proof. By Theorem 4.13, if S then S for fc = 1,..., TV. Moreover, 
for all V e A/ 3 (s,.g,,g,+i)(iP“G) = using the fact that ^ and 

/r' are x*—related, we have, for ^ € T*{V°‘G), 

(/3(7(g,.g,+i)),C) = (/5'(7(g'.g'^j),^x(0) 

(®(7(gfc+i,gfc+2))’ 

and therefore = a{7{gk+i,gk+.)) for fc = 1,..., TV - 2. That is, 7 (gi.g 2 )> 

..., 7 (g„_i,g„) satisfy the discrete second-order dynamics for the discrete Lagrangian 
Ld : Gs ^ M. □ 

Finally, we introduce the notion of Noether symmetry and constants of motion 
for discrete second-order Lagrangian systems and we prove that for all Noether 
symmetry of the discrete second-order Lagrangian system determined by : G 2 —t 
K there is a corresponding constant of motion which is preserved by the discrete 
second-order dynamics. This is a natural extension of discrete Noether symmetry 
for first order systems introduced in [41] and [43]. 

Definition 4.14. A section Z G r(A('P“G')) is said to be a Noether symmetry of 
the discrete second-order Lagrangian system determined by : G 2 —t K if there 
exists a function / G C°°{G) such that 

(d(/r), Z(a“(0)) + /(a“(0) = (0)) + /{^(0) 

for all ^ G Ei, with f = TT-pcQ{p), where ir-pcQ : T*{V°‘G) —>■ V°‘G is the cotan¬ 
gent bundle projection and a“,/3“ : V°‘G —>■ G are the source and target maps, 
respectively, of the Lie groupoid V°'G over G. 

When / = 0, we say is invariant with respect to Z, and the conserved quantity 
is 

F^(m) = (FL±(m),Z), 

and when / 7 ^ 0, we say Ld is quasi-invariant with respect to Z. 

Theorem 4.15. If Z G r{A{'P°‘G)) is a Noether symmetry of a discrete second- 
order Lagrangian system determined by the discrete second-order Lagrangian Ld : 
G 2 —t K, the function Fz : E^ —>■ M given by 

Fziyi) = (d(/x), Z(a“(0)) + /(«“(?)) = (/3 (m), (0)) + /(/3“(e)), 

is a eonstant of motion where ^ = n-pacih')- That is, */ 7 (gi,g 2 )) • ■ • )7(giv-i.gN) ^ 
T*{V°‘G) satisfy the diserete second-order dynamics then, 

-F!z(7(gfe ,gfe + i)) -Fz (7(gfc+i ,gfe+ 2 )) 


for k = 1,...,N -2. 

Proof. If 7(gi,g2)) • ■ • )7(gjv-i.giv) satisfy the discrete second-order dynamics, then 
^(7(g,,g.+i)) =fi(7(gfc+i.g,+2)) where 7(g,.g,+i) G Ep fov k = 1,..., N- 1. Therefore 
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Pzin{gk,gk+i)) 9k,9k+i )),-2'(^“(fffe,5fc,5fc+i))) +/(/3“(Sfe,5fe,5fc+i)) 

= 0h{gk,gk+i))^ Zigk+i)) + /(fffc+i) 

{^{l(gk+i,gk+ 2 ))’ (^fc+l) 5 fe+l) 5/0+2))) 

+ /(a“(5fc+i,5fc+i,5fc+2)) 

Fz {"/{gk+i ,gfc+2)) 


for fc = 1,..., iV — 2. □ 

4.5. Discrete second-order constrained mechanics. Consider a discrete second- 
order constrained systems, that is, given a Lie groupoid C? =4 Q one consider the 
discrete second-order Lagrangian : iV —>■ K defined on a submanifold N of the 
set of composable elements 6 * 2 ,. The submanifold N implies that the dynamics is 
restricted. 

From Theorem 4.3 


is a Lagrangian submanifold where iN : N ^ 'P°‘G is the inclusion from N to 'P°‘G. 

The Lagrangian submanifold, is an affine bundle over N taking the pro¬ 
jection TT-pcQ Ff, the restriction of the cotangent bundle projection 

TT'pcQ : T*{V°'G) —>■ V^G to this Lagrangian submanifold. 

Suppose that the constraint submanifold N C G 2 is given by 

N = {( 51 , 52 ) e G 2 I $'^( 51 , 52 ) = 0, with A e /}, 

where is a family of real functions defined in a neighborhood of N and I 

is an index set. Then, an element fj, S with ( 51 , 52 ) = Ie« (m), can be 

written as 

5 = dLd{gi, g2) + g2) = d{Ld + XA^^)igi,g2) € 

where : ^2 —>■ K is an arbitrary extension of : N ^ M.. 

In this sense, can be locally seen as the space consisting of the elements 
( 51 , 52 ) G N together the Lagrange multipliers constraining ( 51 , 52 ) to N. 

Therefore, by Theorem 4.5, the sequence ( 51 , 52 , ■ ■ •, 5 Af, A^,..., A^“^), is a 
solution of the discrete second-order constrained Lagrangian system determined by 
Ld = Ld + AA<i>^ with ( 5 j, 5 j+i) £ N for j = 1,..., iV — 1 if it satisfies 

0 =(^[gk+i), DiLd{gk+i, gk+ 2 ) + Di^d {gk+i, gk+ 2 ) + D 2 Ldigk,gk+i) 

+ X\D 2 ^d{gkj gk+i)^ ~ (^{gk),DiLd{gk,gk+i) + [gk^ gk+i) 

+ D2Ld{gk-i,gk) + x\ ^D2^d{gk-i^gk)^, 

for fc = 2 ,..., iV — 2 and X a vector field on G. 
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Therefore, the sequence (( 71 , 52 , • ■ ■ ,9n, A^, ..., satisfies 

0 =‘^d(.9k,9k+i), for all A e/, fc = 1,..., - 1; 

0 (^DiLd{9k, 9k+i) + X'XDi^^{gk, 9k+i) + D2Ld{9k-i, 9k) 

+ ^D2^di9k-i,9k)'^ + {fgk+i ° i)* (^DiLd{9k+i, 9k+2) 

+ X'^^Di^^{gk+i, 9k+2) + D2Ld{gk, 9k+i) + ^A^'^^di9k, 9k+i)^, 
for fc = 2,..., iV - 2. 

Remark 9. When the Lie groupoid is a Lie group, we obtain the second-order 
Euler-Poincare equations for systems with constraints (see [17] for example) 

0 =^di9k,9k+i), 0 = ^di9k-i,9k), 0 = $;^( 5 fc+i, 5 fc+ 2 ) for all A e/; 

(^^i^d{9k,9k+i) + X'XDi^^{gk,gk+i) + D2Ld{gk-i,9k) 

+ x’^^D 2 ^d i9k-i, 9k)^ - (^DiLd{9k+i,9k+2) + X'^^Di^-^{gk+i, gk+ 2 ) 

+ D2Ld{gk, 9k+i) + X\D2^di9k, 9k+i)^, 
for fc = 2,..., TV - 2. 

When the Lie groupoid is the Banal groupoid, we have 

0 ('?/c,5fc+i,g/c+2), 0 = $;^(gfc_i,5fc,5fc_^i), 0 ^ <^d{9k-2,qk-i,qk) 

0 =DiLd{qk, qk+i, qk+ 2 ) + ( 9 fe, 5 fc+i, 5 fc+ 2 ) + D 2 Ld{qk-i,qk, qk+i) 

+ x\ ^D2^diqk-i,qki qk+i) + D3Ld{qk-2,qk-i, qk) 

+ A^^^L's^’d (9/c-2,5fe-i,5fc), 

for all A G J and for k = 2,..., N — 2. The equations given above are the discrete 
second-order Euler-Lagrange equations for systems with second-order constraints 
(see [14] for example). o 

5. Conclusions and future research. In this paper, we have developed a gen¬ 
eralized theory of discrete second-order Lagrangian mechanics from a variational 
point of view and we have shown how to apply this theory to the construction of 
variational integrators for some interesting examples of optimal control problems 
of mechanical systems. After that, we have shown how Lagrangian submanifolds 
of a symplectic groupoid (cotangent groupoid) give rise an intrinsic way to discrete 
dynamical second-order systems, and we have studied the geometric properties of 
these systems from the perspective of symplectic and Poisson geometry. Finally, 
we have developed the reduction by Noether symmetries, and we have studied the 
relationship between the dynamics and variational principles for these second-order 
variational problems. 

In [3] we have studied optimal control problems for nonholonomic mechanical 
systems as second-order constrained variational problems. Let I? be a non-integrable 
distribution defined by the nonholonomic constraints of some mechanical systems. 
We define the submanifold of TV by 

:= {u G TV \ V = 7 ( 0 ) where 7 : / —>■ I? is admissible}. 


(51) 


34 


LEONARDO COLOMBO AND DAVID MARTIN DE DIEGO 


and where we choose coordinates on and where the inclusion on 

TV, ix)( 2 ) : V^'^^ ^ TV is given by 

^DO)(9^2/^r) = (<?*, y“,pL(9)y“,y“)- 

Therefore, is locally described by the constraints on TV 

= 0 . 

The optimal control problem is determined by a smooth function L : V^^'> —?> K 
given a cost functional as in Section 3. To derive the equations of motion for L 
one can use standard variational calculus for systems with constraints defining the 
extended Lagrangian L, 

L = L + \M-f^o.y")- 

In a future work we would like to build variational integrators as an alternative 
way to construct integration schemes for the type of optimal control problems stud¬ 
ied in [3]. Since the space is a subset of TV we can discretize the tangent 
bundle TV by the cartesian product V x V. Therefore, our discrete variational 
approach for optimal control problems of nonholonomic mechanical systems will 

be determined by the construction of a discrete Lagrangian Ld ■ V^ —)■ K where 
( 2 ) 

Vj^ is the subset ofVxV locally determined by imposing the discretization of the 
constraint g* = pl^{q)y°‘. For instance we can consider 


V 


( 2 ) _ 
d 


(9o.2/o>9i>2/“) ^'DxV 




Now the system is in a form appropriate for the application of discrete variational 
methods for constrained systems developed in this work from both, variational and 
geometrical points of view as in sections 3.2 and 4.5. 


Appendix A: Higher-order tangent bundles. In this Appendix we recall some 
basic facts of the geometry of tangent bundle theory. For more details see [21, 37]. 

Let Q be a differentiable manifold of dimension n. It is possible to introduce an 
equivalence relation in the set (7^(11, Q) of fc-differentiable curves from K to Q. By 
definition, two given curves in Q, 71 (t) and 72 (t), where t G {—a,a) with a G K 
have contact of order k at qo = 71 (0) = 72 ( 0 ) if there is a local chart {(p,U) of Q 
such that qo G U and 

for all s = 0,...,k. This is a well defined equivalence relation in Q) and the 

equivalence class of a curve 7 will be denoted by [ 7 ) 0 ^^ The set of equivalence classes 
will be denoted by T^^'>Q and it is not hard to show that it has a natural structure 
of differentiable manifold. Moreover, Tq : T^^^Q Q where Tq = 7 ( 0 ) is a 

fiber bundle called the tangent bundle of order k of Q. 

From a local chart = (g®) on a neighborhood U oi Q with i = 1,... ,n, it 

is possible to induce local coordinates (g^°\ g*^^\ ..., g^^^) on T<-k)u = (r^ )-\U). 
The standard convention is, g*-®^ = g®, g^^^ = g® and g*-^^ = g®. 


Appendix B: Discrete Mechanics. This appendix briefly reviews some key re¬ 
sults of discrete mechanics (see Marsden and West [51] for more details). 
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A.l. Discrete Lagrangian Mechanics. A discrete Lagrangian is a differentiable 
function : Q x Q —> M, which may be considered as an approximation of the action 
integral defined by a continuous regular Lagrangian L: TQ —)• M. That is, given a 
time step h > 0 small enough, 

hd{qo,qi)- [ L{q{t),q{t)) dt, 

Jo 

where q{t) is the unique solution of the Euler-Lagrange equations for L with bound¬ 
ary conditions g( 0 ) = qo and q{h) = qi. 

We construct the grid {t/, = kh \ k = 0,...,N}, with Nh = T and define 
the discrete path space Vd{Q) ■= {qd ■ {tk}k=o *3}- identify a discrete 
trajectory qd G VdiQ) with its image qd = {9fe}^o '''^here qk := qd{tk)- The discrete 
action Ad ■ 'PdiQ) —t K along this sequence is calculated by summing the discrete 
Lagrangian on each adjacent pair and defined by 

AT-l 

Ad{qd) = Ad{qo,-;qN) ■■= '^hd{qk,qk+i)- (52) 

fc =0 

We would like to point out that the discrete path space is isomorphic to the smooth 
product manifold which consists on A -|- 1 copies of Q, the discrete action inherits 
the smoothness of the discrete Lagrangian and the tangent space Tq^VdiQ) Qd is 
the set of maps Vq^ : —>■ TQ such that tq o Vqj^ = qd which will be denoted 

by Vq^ = {iqk,Vk)}^=o, where tq : TQ Q is the canonical projection. 

For any product manifold Q1XQ2, 'r^qj^q^^{QixQ2) ~ Qi xT*^(52, for qi G Qi 
and q 2 G Q 2 where T*Q denotes the cotangent bundle of a differentiable manifold 
Q. Therefore, any covector a G 52 ) (Qi X Q 2 ) admits an unique decomposition 
a = ai -I- a2 where at G T*_Qi, for i = 1 , 2 . Thus, given a discrete Lagrangian 
we have the following decomposition 

d]Lj;((7o, qi) = T>ihd{qo, qi) + D 2 hd{qo, qi), 

where Di'Ld{qo,qi) G T*^Q and D 2 Ld(qo, 9 i) G T*^Q. 

The discrete variational principle, or Cadzow’s principle [10], states that the 
solutions of the discrete system determined by must extremize the action sum 
given fixed points qo and q^- Extremizing Ad over qk with 1 < fc < A — 1, we obtain 
the following system of difference equations 

Dihdiqk, qk+i) + D 2 'Ldiqk-i,qk) = 0. (53) 

These equations are usually called discrete Euler-Lagrange equations. Given a 
solution {ql,}kez of eq.(53) and assuming the regularity hypothesis (the matrix 
{Di2hd{qk, Qk+i)) is regular), it is possible to define implicitly a (local) discrete 
flow : Uk C Q X Q ^ Q X Q, hy i:i^^{qk-i,qk) = {qk,qk+i) from (53) where Uk 
is a neighborhood of the point {ql_i,qk)- 

Let us define the discrete Lagrangian 1-forms 0^^ : Q x Q ^ T* {Q x Q) hy 

: {qk,qk+i) ' — D 2 hdiqk,qk-i-i) dqk+i, (54a) 

: {qk,qk+i) '— -Di'Ld{qk,qk-i-i) dqk. (54b) 

Then, the discrete flow preserves the discrete Lagrangian form 

^^Ld(9/c,gfc-i-i) = -d0]Ld = “d0L_j = Ilii:) 2 Ld(gfc,<?fc-i-i) dqk A dqk+i- 


(55) 
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Specifically, we have 

B.2. Discrete Hamiltonian Mechanics. Introduce the right and left discrete 
Legendre transformations FL^ : Q x Q ^ T*Q hy 

FL;!" : {qk,qk+i) ' — (g^+i, £> 2 Ld(gfc, gfe+i)), (56a) 

FL“ : {qk,qk+i) '—5^ (^fc,-DiLd(gfe,Qfc+i)), (56b) 

respectively. Then we find that the Eq. (54) and (55) are pull-backs by these 
maps of the Liouville 1-form Aq and the canonical symplectic 2-form uiq, on T*Q, 
respectively, as follows: 

0± = (FL±)*Aq, = (FL±)*a.Q. 

Let us define the momenta 

Pk,k+i = -Di'LdiQk, Qk+i), Pk,k+i = ^2Ld(gfc, (?fc+i). 

Then, the discrete Euler-Lagrange equations become simply Pk-i k — Pk k+i- 
defining 

Pk = Pk-l,k — Pk,k+l^ 

one can rewrite the discrete Euler-Lagrange equations as follows: 

Pk = -Dihdiqk, Qk+i), 

Pk+i = D2Ld(gfc,gfc+i). 

Furthermore, define the discrete Hamiltonian map : T* Q ^ T* Q hy 

■ {qk,Pk){qk+i,Pk+i)- (58) 

Then, one may relate this map with the discrete Legendre transforms in Eq. (56) 
as follows: 

EL, =FL+o(FL-)-b (59) 

Furthermore, one can also show that this map is symplectic, i.e., 

(EL,)*a;Q = ujQ. 

This corresponds to the Hamiltonian description of the dynamics defined by the 
discrete Euler-Lagrange equations introduced by Marsden and West in [51]. No¬ 
tice, however, that no discrete analogue of Hamilton’s equations is introduced here, 
although the flow is now on the cotangent bundle T*Q. 


Appendix C: Prolongation of Lie algebroids and Mechanics on Lie al- 
gebroids. In this Appendix we recall the definition of the prolongation of a Lie 
agebroid ta ■ A ^ Q over its projection map and the Euler-Lagrange equations on 
Lie algebroids. Further details can be found in [38], [44] and [46]. 
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C.l. Prolongation of a Lie algebroid. Let (A, I-,-],/;)) be a Lie algebroid of 
rank n over Q with projection ta ■ A ^ Q. 

The prolongation of A over its canonical projection, also called A tangent bundle 
to A, is defined to be 

T^^A = U GAxTaA\ p{a') = {TaTA){Va)} 

where Tta : TA ^ TQ is the tangent map to ta- 

In fact, T'^^A is a Lie algebroid of rank 2n over A where : T'^^A —>■ A is the 
vector bundle projection given by T^^\a',Va) = TTA(Pa) = o, and the anchor map 
is Pi := pr 2 : T'^^A TA, the projection over the second factor (see [38] and [46] 
for more details). 

If we now denote by (a, a', Va) an element (o', Va) G T'^^A where a G A and where 
V is tangent, we rewrite the definition of the prolongation of the Lie algebroid as 
the subset oi Ax AxT A given by 

T'^^A = {{a,a' ,Va) G AxAxTA \ p{a') = {TTA){va),Va G TaA and TA{a) = TA{a')}. 
In this sense, if {a,a',Va) G T'^'^A, then pi(a, a', Va) = (a,Va) G TaA, and the 
projection is given by T^\a,a',Va) = a G A. 

The prolongation of A over ta takes the role of TTQ in standard Lagrangian 
mechanics. 

Along the paper when A = AG —>■ Q is a Lie algebroid associated to a Lie 
groupoid we have that T'^^^AG = A{V°‘G) 


Example. Let g be a finite dimensional real Lie algebra, g is a Lie algebroid over a 
single point Q = {(?}. Using that the anchor map of g is zero we obtain that 

G g X Tg} ~ g X g X g ~ 3g. 

The anchor map of T'^‘>g is Pi{^i,^ 2 ,^ 3 ) = (CijCs) G T^^g and the Lie bracket is 
defined by = (K,f?],0). 


C.2. Mechanics on Lie algebroids. (see [44]) 

Let A be a Lie algeroid over Q. We take local coordinates (q*) on Q and a local 
basis {cq,} of sections of the vector bundle ta ■ A ^ Q with a = 1,..., n, then we 
have the corresponding local coordinates on an open subset t^^{U) of A, 

{U is an open subset of Q), where y°‘{a) is the a-th coordinate of a G A in the given 
basis i.e., every a G A is expressed as a = ?/^ei(r^(a)) + ... + y'^en(TA{a)). 

Such coordinates determine local functions on Q which contain the local 

information of the Lie algebroid structure, and accordingly they are called structure 
functions of the Lie algebroid. They are given by 


d 


P(ea)=pL^ and |e„, C/j] = 


These functions should satisfy the relations 

dp 

and 


, -rp _ J dpa^ ^ ^ 

finj Pp dqi 


(60) 


(61) 


E 

cyclic{oi,^,'y) 






(62) 
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which are usually called the structure equations. 

Given a Lagrangian L : A —> K, we fix two points qo, qx in the base manifold 
Q, then we look for admissible curves ^ / C M — > A, (i.e., curves on A such that 
satisfying the variational principle 

0 = 6 [ L{^{t))dt. 

Jo 

The infinitesimal variations are 6 ^ = r]^, for all time-dependent sections rj S 
r(TA), with 77 ( 0 ) = 0 and r]{T) = 0; where is a time-dependent vector field on 
A, the complete lift, locally defined by 

From this variational principle (see [44] for more details) one can derive the 
Euler-Lagrange equations on Lie algebroids 


Jt \^) 




dq^ 

dt 


= PaV 


Appendix D: The Cayley map. The Cayley map cay : g —?► G is defined by 

cay(e)=(e -0 (e + 0 


where e is the identity element of G. The Cayley map is valid for a class of qua¬ 
dratic groups (see [27] for example) that include the most interesting Lie groups in 
mechanics and the one studied in this paper, SO(3). Its right trivialized derivative 
and inverse are defined by 


dcay,^ 2/ = (e - I) ^ + |) ’ ^ ^ ^ “ f) ^ I) ' 


D.l. The Cayley map for iS'0(3). The group of rigid body rotations is repre¬ 
sented by 3 X 3 matrices with orthonormal column vectors corresponding to the 
axes of a right-handed frame attached to the body. On the other hand, the algebra 
so(3) is the set of 3 x 3 antisymmetric matrices. A so(3)-basis can be constructed 
as {ei, 62 , 63 }, Bi € so(3) using the hat map • : —>• so(3), where {ei, 62 , 63 } is the 
standard basis for (see [29] for example). Elements ^ S so(3) can be identified 
with the vector w S through ^ = 00 “ Ba, or f = to. Under such identification the 
Lie bracket coincides with the standard cross product, i.e., ad^ p = a; x p, for some 
p € R^. Using this identification we have 

+iTinTF (“+y) ’ 

where L 3 is the 3x3 identity matrix. The right trivialized derivative and inverse 
are expressed as the 3x3 matrices 

dcay^ = (2-^3 + w), dcay,^ = J 3 - - -^ . (64) 
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